Digital predistortion in varying operating conditions

ABSTRACT

Disclosed are implementations for digital predistortion of signals provided to a radio frequency (RF) transmission path configured to transmit radio signals in a plurality of subbands within a spectral range, including a method that includes configuring a digital predistorter for predistorting signals comprising arbitrary spectral content within the spectral range. The configuring includes acquiring data samples representing operation of the RF transmission path to transmit radio signals in different subbands, each sample including a digital input signal representing spectral content concentrated in a respective subband, and updating parameters of the digital predistorter according to the acquired data samples to mitigate non-linear characteristics of the RF transmission path. The method further includes receiving a further input signal representing spectral content in a particular subband within the spectral range, and using the configured predistorter to process the further input signal to yield a predistorted signal for providing to the RF transmission path.

CROSS-REFERENCE TO RELATED APPLICATIONS

This application is a continuation of U.S. application Ser. No. 16/422,448, entitled “DIGITAL PREDISTORTION IN VARYING OPERATING CONDITIONS” and filed May 24, 2019, which claims the benefit of U.S. Provisional Application Nos. U.S. 62/676,617, entitled “Parameterized Universal DPD Compensators,” U.S. 62/676,677, entitled “Digital Predistortion System with Adaptability to Fast Abrupt Changes of Environmental Characteristics,” and U.S. 62/676,682, entitled “Digital Predistortion System Implementation Using Data Samples Selected According to Variety Metrics,” all filed May 25, 2018, the contents of all of which are incorporated herein by reference in their entireties.

BACKGROUND

The present disclosure relates to approaches to digital predistortion, and more particularly to one or more techniques that alone or in combination are applicable to digital predistortion in quickly varying operating conditions.

Power amplifiers, especially those used to transmit radio frequency communications, generally have nonlinear characteristics. For example, as a power amplifier's output power approaches its maximum rated output, nonlinear distortion of the output occurs. One way of compensating for the nonlinear characteristics of power amplifiers is to ‘predistore’ an input signal (e.g., by adding an ‘inverse distortion’ to the input signal) to negate the nonlinearity of the power amplifier before providing the input signal to the power amplifier. The resulting output of the power amplifier is a linear amplification of the input signal with reduced nonlinear distortion. Digital predistorted devices (comprising, for example, power amplifiers) are relatively inexpensive and power efficient. These properties make digital predistorted power amplifiers attractive for use in telecommunication systems where amplifiers are required to inexpensively, efficiently, and accurately reproduce the signal present at their input. However, the computation effort associated with conventional digital predistortion (DPD) adaptation processes is substantial, which adversely effects the predistortion accuracy and robustness due to the limits the computational efforts may place on the number and speed at which DPD coefficients can be derived.

There are several challenges for DPD adaptation implementations involving bi-directional wireless communication (e.g., between a base station (BS) and a user equipment (UE). For example, a UE typically transmits (uploads data) as needed and with only sufficient output power to save battery power. The adjustment of output power level is achieved by changing, for example, the baseband input level, RF gain, and power supply voltage of a power amplifier. Furthermore, the UE transmitter may be controlled (e.g., at the command of the BS) to make use of varying frequency ranges (e.g., channels) for transmission. Yet further, the mode of operation and parameters, such as power level, power supply voltage level, output load of PA (e.g., VSWR) can rapidly change in UE devices. The combination of such variations in operating conditions can change the characteristics of the transmit chain, and more particularly change the non-linear nature of the transmit chain that is being compensated for by the DPD.

Therefore, there is a need for a DPD approach(es) that can provide high performance DPD (e.g., measured by linearity, error rate, spectral mask violation, etc.) over a substantial range of operating conditions, as well being able to adapt to rapid changes of such operating conditions.

SUMMARY

Disclosed are adaptive digital predistortion implementations that allow the wireless device (generally a UE, but also a base station) to perform quick adaptation (μs scale) of its digital predistortion functionality in response to changes to its operating conditions/points (including changes to environmental conditions). Changes to a wireless device's operating conditions can result due to, for example, 1) the sub-carrier within a component carrier changes frequently (e.g., every time slot, with such time slots having durations in tens of μs) and its spectral content can hop in spectrum rapidly, 2) the UE transmitter needs to adjust average output power by either increasing or reducing the DC supply voltage, 3) changes to other environmental parameters such as temperature and VSWR (such changes may also happen frequently), and/or 4) changes due to the fact that 5G millimeter wave spectrum band is much wider than any existing sub-6 GHz bands, and its spectrum ownership is split between many operators (who operate different bands within that spectral range and may possibly use different technologies). With respect to the latter factor, it is noted that for new millimeter wave spectrum bands, the frequency coverage of a UE device can easily exceed 1 GHz, and as a UE moves from one cell to another, or even as it assigned different cellular resources within a particular cell (by a particular operator), the UE operating points can frequently and abruptly change, thus requiring its digital predistortion functionality to change in accordance with the changes to the operating points.

To address the various problems discussed above, several solutions and approaches to facilitate fast, reliable and robust digital predistortion adaptation techniques are proposed. In one example approach, digital predistortion implementations are provided that are configured for adaptability to fast/abrupt changes in spectral localization, average power, VDD, temperature, VSWR, etc. This adaptability to fast/abrupt changes to the operating points of a wireless device can be achieved, for example, through a DPD system implementation configured to derive an adapted, expanded, set of basis functions, from a first set of DPD basis functions based on operating conditions associated with the wireless device and/or the predistorter system.

In another example approach, a wireless communication device with a universal digital predistortion system is implemented, that is configured to collect training samples from multiple narrow RF bands within a wideband RF spectral range, and optimizes a universal digital predistorter operable over the wideband RF spectral range.

In another example approach to facilitate quick and robust adaptations of digital predistortion functionality, a digital predistortion (DPD) system implements digital predistortion adaptation using data samples selected according to variety metrics. The system may be configured to collect data (e.g., for storage in a variety matrix) distributed over a channel, select a subset of the data that is more uniformly distributed than at least some other subsets of the data, and compute, based on the selected subset of data, digital predistortion parameters to weigh basis functions used to control digital predistortion operations.

Accordingly, in some variations, a method for digital predistortion is provided that includes collecting data for a wireless device comprising a transmit chain that includes at least one power amplifier that produces output with non-linear distortions. The data comprises data records with each data record including input signal data and corresponding output signal data for the transmit chain and is further associated with a set of operating conditions values. Each data record is also associated with a variety measure representative of an extent of distinguishability of the respective data record from other collected data records. The method further includes selecting at least some of the data records according to a variety criterion computed based on the variety measure for the each data record, computing, based on the selected at least some of the data records, digital predistortion parameters used to control a digital predistorter of the wireless device, and applying the digital predistorter, configured with the computed digital predistortion parameters, to subsequent input signals.

Embodiments of the method may include at least some of the features described in the present disclosure, including one or more of the following features.

The set of operating conditions for the each data record may include one or more of, for example, channel information associated with a corresponding channel used by the wireless device to process the respective data record, a corresponding temperature, a corresponding average power, a corresponding power supply voltage (V_(DD)) associated with the at least one power amplifier of the transmit chain, a corresponding Voltage Standing Wave Ratio (VSWR) associated with the wireless device, and/or age of the each data record.

Collecting the data for the wireless device may include collecting the data in one or more spectral channels, with the one or more spectral channels located within a wide band spectral range, and with each of the one or more spectral channels having a respective bandwidth smaller than a range bandwidth of the spectral range.

Each of the data records may include one or more samples obtained by the wireless device under a particular set of operating conditions.

The each of the data records may be maintained in a separate column of a records matrix data structure.

Selecting the at least some of the data records according to the variety criterion computed based on the variety measure for the each data record may include selecting data records associated with respective variety measures indicative of high distinguishability in order to increase diversity of the data records used to compute the digital predistortion parameters.

The method may further include computing, for the each data record, the variety measure based on a distance metric representative of distances between data elements of a particular data record, to respective data elements for at least some of the collected data records.

The method may further include maintaining the computed variety measure for the each data record in a distance vector data structure separate from a data structure to store the collected data records.

The method may further include collecting additional data records, storing a particular one of the collected additional data records in a records data structure in response to a determination that a variety measure computed for the particular one of the collected additional data records is more distinguishable, relative to previously stored data records available at the records data structure, than one of the previously stored data records, and re-computing, at regular or irregular intervals, the digital predistortion parameters based on the stored data samples available in the records data structure, including the stored particular one of the collected additional data records.

Storing the particular one of the collected additional data records may include replacing in the records data structures the one of the previously stored data records with the particular one of the collected additional data records.

The method may further include computing a corresponding variety measure for the particular one of the collected additional data records, including determining for each pair of records comprising the particular one of the collected additional data records and at least one of the previously stored data records a maximal absolute value of difference between the corresponding elements of the data records.

Applying the digital predistorter, configured with the digital predistortion parameters, may include weighing basis functions of the digital predistorter, with the basis functions comprising complex terms applied to values of input samples provided in the data records, and with at least some of the basis functions further comprising real terms applied to values of sets of operating conditions associated with the data records.

The method may further include deriving from a first set of basis functions the basis functions of the digital predistorter, the basis functions of the digital predistorter being an expanded set of basis functions that includes (d+1) groups of basis functions, with d being representative of a number of parameters of the operating conditions, wherein a first group of the expanded set of basis functions includes the first set of basis functions without any modification, and with each of the remaining d groups of the expanded set of basis functions being determined as the product of the first set of basis functions and one of the parameters of the operating conditions.

In some variations a linearization system to perform digital predistortion is provided that includes a transmit chain that includes at least one power amplifier that produces output with non-linear distortions, and a digital predistorter comprising an adaptation section and a compensator. The digital predistorter is configured to collect data comprising data records, with each data record including input signal data and corresponding output signal data for the transmit chain and further being associated with a set of operating conditions values, with the each data record being associated with a variety measure being representative of an extent of distinguishability of the respective data record from other collected data records. The digital predistorter is further configured to select at least some of the data records according to a variety criterion computed based on the variety measure for the each data record, compute, based on the selected at least some of the data records, digital predistortion parameters used to control the digital predistorter, and digitally predistort, according to the computed digital predistortion parameters, subsequent input signals.

Embodiments of the linearization system may include at least some of the features described in the present disclosure, including any of the above method features, as well as one or more of the following features.

The digital predistorter configured to select the at least some of the data records according to the variety criterion computed based on the variety measure for the each data record may be configured to select data records associated with respective variety measures indicative of high distinguishability in order to increase diversity of the data records used to compute the digital predistortion parameters.

The digital predistorter may further be configured to compute, for the each data record, the variety measure based on a distance metric representative of distances between data elements of a particular data record, to respective data elements for at least some of the collected data records.

The digital predistorter may further be configured to collect additional data records, store a particular one of the collected additional data records in a records data structure in response to a determination that a variety measure computed for the particular one of the collected additional data records is more distinguishable, relative to previously stored data records available at the records data structure, than one of the previously stored data records, and re-compute, at regular or irregular intervals, the digital predistortion parameters based on the stored data samples available in the records data structure, including the stored particular one of the collected additional data records.

The digital predistorter configured to digitally predistort the subsequent input signals may be configured to derive from a first set of basis functions an expanded set of basis functions of the digital predistorter, the expanded set of basis functions including (d+1) groups of basis functions, with d being representative of a number of parameters of the operating conditions, with a first group of the expanded set of basis functions including the first set of basis functions without any modification, and with each of the remaining d groups of the expanded set of basis functions being determined as a transformation of the first set of basis functions according to a function of at least one of the parameters of the operating conditions. The digital predistorter may further be configured to weigh the expanded set of basis functions with the computed digital predistortion parameters.

In some variations, another method for digital predistortion of signals provided to a radio frequency (RF) transmission path configured to transmit radio signals in a plurality of subbands within a spectral range, is provided. The method includes configuring a digital predistorter for predistorting signals comprising arbitrary spectral content within the spectral range, the configuring including acquiring one or more data samples representing operation of the RF transmission path to transmit radio signals in different subbands, each data sample including a digital input signal provided to the RF transmission path for transmission and a digital sensed signal representing a sensing of a radio transmission in response to the provided digital signal, the each digital input signal representing spectral content concentrated in a respective single corresponding subband within the spectral range, and updating parameters of the digital predistorter according to the acquired one or more data samples to mitigate non-linear characteristics of the RF transmission path. The method further includes receiving a further input signal for transmission by the RF transmission path, the further input signal representing spectral content in a particular subband within the spectral range, and using the configured predistorter to process the further input signal to yield a predistorted signal for providing to the RF transmission path.

Embodiments of the additional method may include at least some of the features described in the present disclosure, including any of the above features of the first method and the first linearization system, as well as one or more of the following features.

The spectral range, W, may include n RF bands of approximately equal-sized bandwidths.

Acquiring the one or more data samples may include controllably adjusting a wireless device comprising the RF transmission path to operate in one or more of the different subbands within the spectral range, and at each of the one or more of the different subbands, measuring at least one set of samples.

Acquiring the one or more data samples may include measuring during runtime, subsequent to an initial computation of the parameters of the digital predistorter, additional sets of samples at additional subbands within the spectral range, and re-updating the parameters of the digital predistorter further based on at least some of the additional sets of samples.

The digital predistorter may be implemented in a user equipment (UE), and the plurality of subbands may be successively assigned bands for transmission.

The updated parameters of the digital predistorter may be configured to weigh basis functions used to predistort the further input signal, with the basis functions including complex terms that depend on a τ-dimensional complex vector q corresponding to complex input samples to the digital distorter, and real terms that depend on a d-dimensional real vector p representative of values of operating conditions associated with operation of the RF transmission path, with d being representative of a number of parameters of the operating conditions associated with the operation of the RF transmission path.

A total number of the basis functions used to digitally predistort input signals may be proportional to the product of a number of terms, b, in a first set of basis functions, and d+1, such that the total number of basis functions is proportional to b·(d+1).

Acquiring the one or more data samples representing the operation of the RF transmission path in the different subbands may further include collecting for the one or more data samples corresponding operating conditions parameters, with the corresponding operating conditions parameters including one or more of, for example, a corresponding temperature of the RF transmission path, a corresponding average power, a corresponding power supply voltage (V_(DD)) associated with at least one power amplifier of the RF transmission path, a corresponding Voltage Standing Wave Ratio (VSWR) associated with the RF transmission path, age of the collected each of the sets of samples, and/or corresponding channel information in which the one or more data samples is collected.

Using the configured predistorter to process the further input signal may include predistorting the further input signal according to:

${{v\lbrack t\rbrack} = {{u\lbrack t\rbrack} + {\sum\limits_{k = 1}^{n}{x_{k}{B_{k}\left( {{q_{u}\lbrack t\rbrack},{p\lbrack t\rbrack}} \right)}}}}},{{{where}\mspace{14mu}{q_{u}\lbrack t\rbrack}} = \begin{bmatrix} {u\lbrack t + l - 1\rbrack} \\ \vdots \\ {u\lbrack t + l - \tau\rbrack} \end{bmatrix}}$ where v[t] is an output of the digital predistorter at time t, u[t] is an input sample to the digital predistorter at the time t, B_(k) are the basis functions for k=1 to n, q_(u)[t] are stacks of input samples, p[t] are d-dimensional real vectors of operating conditions at the time t, and x_(k) are the computed digital predistortion parameters to weigh the basis functions.

Updating the parameters of the digital predistorter may include deriving, according to an optimization process that is based, at least in part, on the acquired one or more data samples measured in the different subbands, the digital predistortion parameters, x_(k), to weigh the basis functions B_(k).

Deriving according to the optimization process the parameters, x_(k), of the digital predistorter may include adaptively adjusting the digital predistortion parameters, x_(k), through minimization of a regularized mean square error function E(x) according to:

${{E(x)} = {{\frac{\rho}{n}{x}^{2}} + {\frac{1}{N}{\sum\limits_{t \in T}{{{v\lbrack t\rbrack} - {y\lbrack t\rbrack} - {\sum\limits_{k = 1}^{n}{x_{k}{B_{k}\left( {{q_{u}\lbrack t\rbrack},{p\lbrack t\rbrack}} \right)}}}}}^{2}}}}},$ where ρ>0 is a regularization factor, v[t] is an output of the digital predistorter at the time t, and y[t] is the output of an observation receiver coupled to an output of the RF transmission path.

Adaptively adjusting the digital predistortion parameters, x_(k), may include deriving correlation matrices G and L, when n<<100, with n being representative of the number of basis functions, according to:

${G = {\frac{1}{N}{\sum\limits_{t \in T}{{B\lbrack t\rbrack}{B\lbrack t\rbrack}^{\prime}}}}},{L = {\frac{1}{N}{\sum\limits_{t \in T}{\left( {{v\lbrack t\rbrack} - {y\lbrack t\rbrack}} \right)^{\prime}{B\lbrack t\rbrack}}}}},{where}$ ${{B\lbrack t\rbrack} = \begin{bmatrix} {B_{1}\left( {{q_{y}\lbrack t\rbrack},{p\lbrack t\rbrack}} \right)} \\ \vdots \\ {B_{n}\left( {{q_{y}\lbrack t\rbrack},{p\lbrack t\rbrack}} \right)} \end{bmatrix}};$ and deriving the digital predistortion parameters, x_(k), according to:

$\overset{¯}{x} = {\left( {{\left( \frac{p}{n} \right)I} + G} \right)^{- 1}{L.}}$

Adaptively adjusting the digital predistortion parameters, x_(k), may include approximating the digital predistortion parameters, x_(k), according to a stochastic gradient process.

Acquiring the one or more data samples representing operation of the RF transmission path to transmit radio signals in the different subbands may include selecting the one or more data samples from a plurality of data samples according to a variety criterion computed based on the variety measure for the each data record, with the plurality of records being arranged in a variety database structure in which each record corresponds to one of the one or more data samples.

In some variations, an additional linearization system to perform digital predistortion of signals provided to a radio frequency (RF) transmission path configured to transmit radio signals in a plurality of subbands within a spectral range, is provided. The system includes the RF transmission path comprising at least one power amplifier that produces output with non-linear distortions, and a digital predistorter comprising an adaptation section and a compensator. The digital predistorter is configured to configure the compensator to predistort signals comprising arbitrary spectral content within the spectral range, including to acquire one or more data samples representing operation of the RF transmission path to transmit radio signals in different subbands, each data sample including a digital input signal provided to the RF transmission path for transmission and a digital sensed signal representing a sensing of a radio transmission in response to the provided digital signal, the each digital input signal representing spectral content concentrated in a respective single corresponding subband within the spectral range. Configuring the compensator also includes updating parameters of the digital predistorter according to the acquired one or more data samples to mitigate non-linear characteristics of the RF transmission path. The digital predistorter is additionally configured to receive a further input signal for transmission by the RF transmission path, the further input signal representing spectral content in a particular subband within the spectral range, and use the configured compensator to process the further input signal to yield a predistorted signal for providing to the RF transmission path.

Embodiments of the additional linearization system may include at least some of the features described in the present disclosure, including any of the above features for the various methods and for the first linearization system.

The digital predistorter configured to acquire the one or more data samples may be configured to controllably adjust a wireless device comprising the RF transmission path to operate in one or more of the different subbands within the spectral range, and at each of the one or more of the different subbands, measure at least one set of samples.

The digital predistorter configured to acquire the one or more data samples may be configured to measure during runtime, subsequent to an initial computation of the parameters of the digital predistorter, additional sets of samples at additional subbands within the spectral range, and re-update the parameters of the digital predistorter further based on at least some of the additional sets of samples.

The digital predistorter configured to use the compensator may be configured to derive from a first set of basis functions an expanded set of basis functions of the digital predistorter, the expanded set of basis functions including (d+1) groups of basis functions, with d being representative of a number of parameters of operating conditions, with a first group of the expanded set of basis functions including the first set of basis functions without any modification, and with each of the remaining d groups of the expanded set of basis functions being determined as a transformation of the first set of basis functions according to a function of at least one of the parameters of the operating conditions. The digital predistorter may further be configured to weigh the expanded set of basis functions with the updated digital predistortion parameters, to process the further input signal.

The digital predistorter configured to acquire the one or more data samples representing operation of the RF transmission path to transmit radio signals in the different subbands may be configured to select the one or more data samples from a plurality of data samples according to a variety criterion computed based on the variety measure for the each data record, wherein the plurality of records are arranged in a variety database structure in which each record corresponds to one of the one or more data samples.

In some variations, a further method for digital predistortion is provided, that includes obtaining a first set of digital predistortion (DPD) basis functions for controlling operation of a digital predistorter of a wireless device operating on a received at least one input signal directed to a power amplification system comprising a transmit chain with at least one power amplifier that produces output with non-linear distortions, and determining an adapted set of DPD basis functions derived from the first set of DPD basis functions, that depend on operating conditions associated with operation of the wireless device, for controlling the operation of the digital predistorter. The method also includes configuring the digital predistorter with DPD coefficients determined for the adapted set of DPD basis functions based, at least in part, on observed samples of the transmit chain responsive to the at least one input signal.

Embodiments of the further method may include at least some of the features described in the present disclosure, including any of the above features of the first and second methods, and the first and second linearization systems, as well as one or more of the following features.

The basis functions may include complex terms that depend on a τ-dimensional complex vector q corresponding to complex input samples to the digital distorter, and real terms that depend on a d-dimensional real vector p representative of values of the operating conditions associated with the operation of the wireless device on the received at least one input signal, with d being representative of a number of parameters of the operating conditions associated with the operation of the wireless device.

A total number of basis functions used to digitally predistort the at least one input signal may be proportional to the product of a number of terms, b, in the first set of basis functions, and d+1, such that the total number of basis functions is proportional to b·(d+1).

The basis functions may include (d+1) groups of basis functions, with a first group of basis functions including the first set of basis functions without any modification, and with each of the remaining d groups of basis functions determined as the product of the first set of basis functions and one of the operating conditions associated with the wireless device.

Configuring the predistorter with the DPD coefficients may include computing the DPD coefficients according to an optimization process that yields the DPD coefficients, x_(k), as weights to the basis functions that minimize an error function based on the observed samples of the transmit chain and samples of the at least one input signal.

Computing, according to the optimization process, the DPD coefficients, x_(k), may include adaptively adjusting the DPD coefficients, x_(k), through minimization of a regularized mean square error function E(x) according to:

${{E(x)} = {{\frac{\rho}{n}{x}^{2}} + {\frac{1}{N}{\sum\limits_{t \in T}{{{v\lbrack t\rbrack} - {y\lbrack t\rbrack} - {\sum\limits_{k = 1}^{n}{x_{k}{B_{k}\left( {{q_{y}\lbrack t\rbrack},{p\lbrack t\rbrack}} \right)}}}}}^{2}}}}},{where}$ ${q_{y}\lbrack t\rbrack} = \begin{bmatrix} {y\lbrack t + l - 1\rbrack} \\ \vdots \\ {y\lbrack t + l - \tau\rbrack} \end{bmatrix}$ where ρ>0 is a regularization factor, B_(k) are the basis functions for k=1 to n, q_(y) [t] are stacks of input samples at time t, p[t] are d-dimensional real vectors of the operating conditions at the time t, v[t] is an output of a digital predistortion (DPD) compensator at the time t, and y[t] is the output of an observation receiver coupled to an output of the transmit chain.

Adaptively adjusting the DPD coefficients, x_(k), may include deriving correlation matrices G and L, when n<<100, with n being representative of the number of basis functions, according to:

${G = {\frac{1}{N}{\sum\limits_{t \in T}{{B\lbrack t\rbrack}{B\lbrack t\rbrack}^{\prime}}}}},{L = {\frac{1}{N}{\sum\limits_{t \in T}{\left( {{v\lbrack t\rbrack} - {y\lbrack t\rbrack}} \right)^{\prime}{B\lbrack t\rbrack}}}}},{where}$ ${{B\lbrack t\rbrack} = \begin{bmatrix} {B_{1}\left( {{q_{y}\lbrack t\rbrack},{p\lbrack t\rbrack}} \right)} \\ \vdots \\ {B_{n}\left( {{q_{y}\lbrack t\rbrack},{p\lbrack t\rbrack}} \right)} \end{bmatrix}};$ and deriving the DPD coefficients, x_(k), according to:

$\overset{\_}{x} = {\left( {{\left( \frac{p}{n} \right)I} + G} \right)^{- 1}{L.}}$

Adaptively adjusting the DPD coefficients, x_(k), may include approximating the DPD coefficients, x_(k), according to a stochastic gradient process.

The method may further include collecting data records representative of samples of the at least one input signal and respective one or more observed samples of the transmit chain responsive to the samples of the at least one input signal.

Collecting the data records may include collecting the data records in one or more spectral channels, the one or more spectral channels located within a wide band spectral range, with each of the one or more spectral channels having a respective bandwidth smaller than a range bandwidth of the spectral range.

In some variations, a further linearization system to perform digital predistortion is provided that includes a transmit chain comprising at least one power amplifier that produces output with non-linear distortions, and a digital predistorter comprising an adaptation section and a compensator. The digital predistorter is configured to obtain a first set of digital predistortion (DPD) basis functions for controlling operation of the digital predistorter of a wireless device operating on a received at least one input signal, and determine an adapted set of DPD basis functions, derived from the first set of DPD basis functions, that depends on operating conditions associated with operation of the wireless device, for controlling the operation of the digital predistorter. The digital predistorter is also configured to configure the compensator with DPD coefficients determined for the adapted set of DPD basis functions based, at least in part, on observed samples of the transmit chain responsive to the at least one input signal.

Embodiments of the further linearization system may include at least some of the features described in the present disclosure, including any of the above features for the various methods, and for the various linearization systems, as well as one or more of the following features.

The basis functions may include (d+1) groups of basis functions, with a first group of basis functions including the first set of basis functions without any modification, and with each of the remaining d groups of basis functions determined as a transformation of the first set of basis functions according to a function of at least one of the parameters of the operating conditions.

The each of the remaining d groups of basis functions may be determined as the product of the first set of basis functions and one of the operating conditions associated with the wireless device.

The digital predistorter configured to configure the compensator with the DPD coefficients may be configured to compute the DPD coefficients according to an optimization process that yields the DPD coefficients, x_(k), as weights to the basis functions that minimize an error function based on the observed samples of the transmit chain and samples of the at least one input signal.

The digital predistorter may further be configured to collect data records, representative of samples of the at least one input signal and respective one or more observed samples of the transmit chain responsive to the samples of the at least one input signal, in one or more spectral channels, the one or more spectral channels located within a wide band spectral range, with each of the one or more spectral channels having a respective bandwidth smaller than a range bandwidth of the spectral range.

In some variation, yet another method for digital predistortion is provided that includes collecting data for a wireless device comprising a transmit chain that includes at least one power amplifier that produces output with non-linear distortions. The data includes data records with each data record including input signal data and corresponding output signal data for the transmit chain, and further being associated with a set of operating conditions, with the each data record is associated with a variety measure representative of an extent of distinguishability of the respective data record from other collected data records, and with collecting the data including collecting the data in one or more spectral channels, the one or more spectral channels being located within a wide band spectral range, with each of the one or more spectral channels having a respective bandwidth smaller than a range bandwidth of the spectral range. The method also includes selecting at least some of the data records according to a variety criterion computed based on the variety measure for the each data record, and computing, based on the selected at least some of the data records, digital predistortion parameters used to control a digital predistorter of the wireless device, with computing the digital predistortion parameters including obtaining a first set of digital predistortion (DPD) basis functions for controlling operations of the digital predistorter, determining an adapted set of DPD basis functions, derived from the first set of DPD basis functions, that depend on operating conditions values corresponding to the selected at least some of the data records, and configuring the digital predistorter with the digital predistortion parameters determined for the adapted set of DPD basis functions based, at least in part, on the selected at least some of the data records. The method also includes applying the digital predistorter, configured with the computed digital predistortion parameters, to subsequent input signals.

Embodiments of the yet other method may include at least some of the features described in the present disclosure, including any of the above features for the various methods, and/or for the various linearization systems.

In some variations, systems are provided that are configured to perform one or more of the method steps provided above.

In some variations, a design structure is provided that is encoded on a non-transitory machine-readable medium, with the design structure including elements that, when processed in a computer-aided design system, generate a machine-executable representation of one or more of the systems, digital predistorters, transmit chains, and/or any of their respective modules, as described herein.

In some variations, an integrated circuit definition dataset is provided that, when processed in an integrated circuit manufacturing system, configures the integrated circuit manufacturing system to manufacture one or more of the systems, digital predistorters, transmit chains, and/or any of their respective modules, as described herein.

In some variations, a non-transitory computer readable media is provided that is programmed with a set of computer instructions executable on a processor that, when executed, cause the operations comprising the various method steps described above.

Other features and advantages of the invention are apparent from the following description, and from the claims.

BRIEF DESCRIPTION OF THE DRAWINGS

These and other aspects will now be described in detail with reference to the following drawings.

FIG. 1 is a block diagram of a linearization system configured to implement various DPD adaptation approaches as described herein.

FIG. 2 is a diagram illustrating collection of samples from narrow subbands of a wideband spectral range.

FIG. 3 includes graphs comparing the performance results of DPD compensation systems implemented using two different approaches (namely, a wideband sample collection approach, and narrow subbands sample collection approach).

FIG. 4 is a flowchart of an example procedure for digital predistortion for a radio frequency (RF) transmission path configured to transmit radio signals in a plurality of subbands within a spectral range.

FIG. 5 is a flowchart of an example procedure for digital predistortion using adaptive basis functions.

FIG. 6A includes graphs illustrating simulation results for an optimization process performed without ensuring a good variety of samples.

FIG. 6B includes graphs showing simulation results when the optimization process used samples that were subject to a variety-based distribution (randomization) approach to promote sample diversity.

FIG. 7 is a flowchart of an example procedure for digital predistortion implemented using a DPD adaptation process based on a variety-based approach.

Like reference symbols in the various drawings indicate like elements.

DESCRIPTION

Very generally, the following includes a description of a number of techniques that may be used individually, or preferably in combination, to provide a digital predistorter that provides high performance over a wide range of operating conditions, and is able to react to rapid changes between such operating conditions.

These techniques are particularly applicable to radio frequency (RF) user equipment (UE), such as a cellular telephone, that is controlled (e.g., by a base station (BS)) to make use of rapid variations, for example, in spectral positions (e.g., transmission channels or subcarriers), power levels, or output loading. However, it should be understood that the techniques are not limited to implementation on user equipment, and that the techniques may be used on base stations as well. Furthermore, although described in the context of predistortion for over-the-air RF applications, the techniques may be applied to other domains including coaxial cable communication (e.g., for use in DOCSIS standard modems), audio communication (e.g., for linearizing audio, ultrasound, ultra-low frequency), and optical communication (e.g., to mitigate non-linear characteristics of optical communication paths).

A first of these techniques is motivated in part by the varying spectral position of a transmit spectrum. In this technique, generally, a predistorter is configured to process input signals that span a relatively narrow subband within a wider available bandwidth, and is “universal” in the sense that it is largely or completely agnostic to the location and/or width of the particular subband in which the transmitter is operating at any given time. It should be recognized that one challenge for this technique is that data available to adapt the predistorter is sparse, in the sense that feedback available to update the parameters of the predistorter only sense characteristics over a small fraction of the overall available bandwidth at any one time.

A second of these techniques explicitly accounts for changes in certain aspects of the operating conditions by introducing parameters representing those aspects, and incorporating dependency on those parameters into the model structure of the predistorter. Using this approach, a change in one of these aspects involves changing the value of the corresponding parameter of the predistorter model without having to wait for a slower adaptation procedure to react to the change in the aspect of the operating condition. As an example, the dependence on transmission power may be parameterized, and the predistorter can react to changes in power virtually instantaneously. In another example, the basis functions, which are generally a function of the desired input, that are used to predistort input signals can be further adapted based on changes to operating conditions.

In example implementations of the second technique, a method for digital predistortion is provided that includes obtaining a first set of digital predistortion (DPD) basis functions for controlling operations of a digital predistorter, of a wireless device, operating on at least one input signal directed to a power amplification system comprising a transmit chain with at least one power amplifier that produces output with non-linear distortions. The method additionally includes determining an adapted set of DPD basis functions, derived from the first set of DPD basis functions, that depend on operating conditions (e.g., including environmental parameters) associated with operation of the wireless device, for controlling operation of the digital predistorter, and configuring the digital predistorter with DPD coefficients determined for the adapted set of DPD basis functions based, at least in part, on observed samples of the transmit chain responsive to the at least one input signal.

A third of these techniques addresses the possibility that during operation, the range of operating conditions is not made use of in a representative manner, or follows a time pattern that may result in adaptation of the predistorter to degrade its ability to change to a new operating condition that may be underrepresented in its influence on the configuration of the predistorter. Very generally, the technique makes use of a sampling of operational data based on which the configuration is adapted that maintains an adequate diversity over the possible operating conditions, thereby providing good performance immediately after a change from one configuration to another.

In example implementations of the third type of techniques, digital predistortion adaptation is implemented using data samples selected according to variety metrics. In some such implementations, a method for digital predistortion is provided that includes collecting data for a wireless device comprising a transmit chain that includes at least one power amplifier that produces output with non-linear distortions. The collected data includes data records, with each data record including input signal data and corresponding output signal data for the transmit chain, and further being associated with a set of operating conditions values. Each data record is also associated with a variety metric representative of an extent of distinguishability of the respective each data record from other collected data records. Examples of such a variety metrics include distance metrics between elements of different records (large distances are generally indicative of high distinguishability between such different records), with one example of such distance computation being the maximal absolute value of differences between the corresponding elements of the data records. The method additionally includes selecting at least some of the data records according to a variety criterion computed based on the variety metric for the each data record, computing, based on the selected at least some of the data records, digital predistortion parameters used to control a digital predistorter of the wireless device, and applying the digital predistorter, configured with the computed digital predistortion parameters, to subsequent input signals.

The techniques outlined above are described below in the context of a predistortion system 100 shown in block diagram form in FIG. 1. In FIG. 1, u is an input signal that represents a desired signal to be transmitted via a transmission section F 120. The signal u comprises a sequence . . . , u[t−1], u[t], u[t+1], . . . of complex I/Q samples provided to the system 100 in digital form. As introduced above, this input signal in general occupies a relatively small part of an overall available bandwidth of the system. However, through implementation of an adaptation process that uses samples collected over a large number of diverse bands/subbands within the available bandwidth (e.g., the range of frequencies of a set of available communication channels) of the system (i.e., the implementation referred to as the First Technique, as further discussed in greater detail below), the adaptive DPD compensation system 100 can be realized as a universal DPD compensator configured to can operate at different spectral bands (e.g., within a broad spectral range, W, that includes bands operated by correspond to different commercial cellular operators, different available commercial bands, different cellular technologies, etc.). The techniques are also applicable for use with similar systems described in PCT Application No. PCT/US2019/031714, filed on May 10, 2019, titled “Digital Compensation for a Non-Linear System,” U.S. application Ser. No. 16/004,594, titled “Linearization System,” filed on Jun. 11, 2018, and U.S. application Ser. No. 16/004,713, titled “Multi-Band Linearization System,” filed on Jun. 11, 2018, which are each incorporated herein by reference.

A compensator C 110 is electrically coupled to the transmission section F 120 that represents the signal chain from a digital signal (at baseband, or possibly at in intermediate frequency band) to a radio frequency signal provided to, or emitting from, an antenna, including components such a modulators, power amplifier, analog transmission line, and the like, which can introduce non-linear characteristics into the signal path. The transmission section F 120 (also referred to as transmission chain or transmit chain) implements a digital-to-analog subsystem that includes, for example, a DAC/PA combination. In wireless implementations, the transmission section F 120 also includes a modulator (e.g., implemented using a local oscillator) to allow transmission of a resultant analog signal within some RF band within the spectral range W that the system 100 is configured to operate in. The particular band may be selected and/or adjusted (e.g., during hand-off operations) according to availability of cellular resources (e.g., available RF communication bands).

The compensator (also referred to as the digital predistorter) C 110, is configured to compensate for non-linear characteristics of the transmission section F 120 (such non-linear behavior may be caused by operation of one or more of the components of the transmission section 120, such as the power amplifier, the DCA, the modulator, etc.) In particular the predistorter receives the desired signal u and outputs a corresponding predistorted signal v, which is a function of u, where this function is parameterized by configuration parameters x. A sensing section S 130 (also referred to as the observation path) senses the output of the transmission section F 120, and provides a sensed signal y in digital form. Ideally, the sensed signal y is the same (or a scaled, e.g., amplified, version) as the desired input u if the predistorter 110 is able to exactly invert all the non-linear effects of the transmission section. The sensing section S 130 implements an analog-to-digital sub-system, representing real time output measurement (e.g., a demodulator/ADC combination) of the output, w, of the transmission section F 120. The block S 130 produces observed digital samples denoted as the signal y. In some embodiments, at least a batch of the observed samples is stored in a memory device (e.g., a buffer or a bank of buffers) that is in electrical communication with at least one of the subsystems of the system 100.

An adaptation section A 140 (also referred to as an adapter) receives the desired signal u and the sensed signal y and determines the predistortion parameters x that are used in order to approach the ideal situation in which y is the same as u. As will be discussed in greater detail below, in some example embodiments, the system 100 is configured to derive and apply a predistortion process that uses adaptive basis functions that are based on the operating conditions (e.g., environmental characteristics) of the system 100. In at least some combinations of techniques that explicitly model aspects of the operating conditions, the adaptation section receives operating parameters p, for example provided by a controller of the system, or sensed by sensing devices/sections 150 (schematically depicted in FIG. 1 as sensing characteristics of the transmission section 120, recognizing that this is only one possible source of such operating conditions). Such sensing devices can include, for example, a temperature sensor, current and/or voltage sensors, RF power meter, circuit and/or processor-based devices that derives operating condition parameters (such as the voltage standing wave ratio, or VSWR, values indicative of records' age, such as timestamps), etc.

More particularly, samples of the intermediary signal v and the observed samples y are provided to the adaptation section 140 A, which implements an adaptation process (e.g., an optimization process, a coefficient selection process to select pre-determined coefficients from a look-up table, a coefficient interpolation process, etc.) to derive a vector x of the weighing coefficients of the basis functions applied by the compensator C 110, based on the data obtained for v, and y (and, in some instances, u). The observed samples y, as well as the intermediary samples v and/or the input samples u, may all be stored in a buffer (not shown, which may be part of any of the sub-systems of the system 100). The content of the buffer (which may be arranged as a matrix, or may be arranged according to some other appropriate data structure) may be occasionally updated as new samples, that are deemed to be more useful for the adaptation process (because they are more recent and/or meet some data criterion such as data uniformity) become available.

In at least some implementations, the compensator (predistorter) C 110 implements a non-linear transformation that may be represented as:

${v\lbrack t\rbrack} = {{u\lbrack t\rbrack} + {\sum\limits_{k = 1}^{n}{x_{k}{{B_{k}\left( {{q_{u}\lbrack t\rbrack},{p\lbrack t\rbrack}} \right)}.}}}}$ In this equation, q_(u)[t] is a window of τ past input values

${q_{u}\lbrack t\rbrack} = \begin{bmatrix} {u\left\lbrack {t + l - 1} \right\rbrack} \\ {u\left\lbrack {t + l - 2} \right\rbrack} \\ \vdots \\ {u\left\lbrack {t + l - \tau} \right\rbrack} \end{bmatrix}$ and p[t] is a vector (or scalar) representing aspects of the current operating condition. Note that in versions that do not explicitly model such aspects, the dependence on p[t] may be omitted. The parameters x shown in FIG. 1 comprise a vector of n complex-values parameters x_(k).

The terms B_(k) (q, p) are referred to as “basis functions,” which can output scalar complex valued function of the inputs q and p. A variety of such functions may be used, for example, as described in PCT Application No. PCT/US2019/031714, filed on May 10, 2019, titled “Digital Compensation for a Non-Linear System.” Note that the number of samples in q_(u) depends on the basis functions used according to the maximum lag τ used by any of the basis functions. Very generally, in some embodiments, the basis functions B_(k) may be selected during operation according to the operating conditions, and the adaption section A 140 modifies (i.e., adapts) the parameters x (which weigh the basis functions) as new data samples v and y are made available to it, e.g., in an incremental updating process. One or more such processes are described later in this document. Thus, the adaptive DPD compensation system 100 may be configured to adaptively determine DPD coefficients (e.g., selected from a lookup tables/matrices containing pre-determined DPD coefficients associated with different operating points, or dynamically derived according to current and previously obtained samples).

As noted, in some example approaches, the basis functions may be adaptable functions that vary according to operating conditions, thus providing a larger number and variety of basis functions beyond an initial set of carefully pre-selected basis functions. Thus, in some examples, the adaptability of the basis functions used for the predistortion operations can be made dependent on a vector p representing d operating conditions that is introduced into the predistortion model by essentially introducing a linear dependence on each of these terms. This can be represented by a replacement of

-   -   x_(k)B_(k) (q, p)         with         x _(3(k−1)+1) B _(k0)(q)x _(3(k−1)+2) B _(k0)(q)p ₁ + . . . +x         _(3(k−2)+d+1) B _(k0)(q)p _(d).

Therefore, the number terms in the adapted parameters x goes up by a factor of d+1. In this example, the complexity of basis functions remains unchanged because the same basis function B_(k0) is used for a group of (d+1) terms. In other implementations, different basis functions may be used for each of the aspects of the operating conditions, such that these basis functions may depend, in part, on contemporaneous environmental parameters values measured/determined during the adaptation process, with the number of total basis functions used to digitally predistort the input samples being dependent on the number of operating conditions used. For example, if two environmental parameters (e.g., temperature and VDD) are measured for collected samples, the number of basis functions used for the digital predistortion process will be three (3) times the number of basis functions used in implementation that do not incorporate environmental parameters (or other types of operation conditions) into the adaptation process.

As noted above, and as will be discussed in greater detail below, the system 100 may also be configured for more robust and efficient adaptability through implementation of a variety matrix adaptation approach (the “third technique” referred to above). In such an approach, past samples of (q_(u)[t],u[t]−y[t], p[t]), or, for certain adaptation procedures (q_(y)[t], v[t]−y[t], p[t]), are not used with a uniform distribution over past points in time. One implementation of such non-uniform sampling is to maintain a buffer V of a subset of past samples. This buffer of samples (q^((i)), v^((i)), p^((i))) is maintained by the system to represent a diversity of past operating conditions, without overrepresenting any particular condition. Generally, when a new time sample (q[t],v[t], p[t]) is made available to the adaptation section, the adaption section A 140 determines whether to retain this time sample in the buffer V, and if so, which current entry in the buffer to eject.

One general approach to determining whether to eject a current sample and replace it with the new time sample is to do so if such a replacement can increase the diversity or uniformity of the sampling of operating conditions in the buffer. For example, an easy-to-compute distance between data samples is defined as the maximum absolute difference between components of the data sample, for each data i in the buffer. In general, it is desirable to have data points in the buffer have as large as possible a minimum distance between data points in order to cover a larger hyperspace, thus promoting diversity in the overall buffer. The new time sample similarly has a minimum distance h_(i) to the data points in the buffer. The data sample with the smallest value of min(r_(i), h_(i)) is then ejected and replaced with the new sample. In the adaptation procedure, the updated buffer V is used to update the predistortion parameters x.

As noted, the adaptation section A 140 is configured to compute the DPD coefficients that weigh basis functions applied to the input signal u[t] by the compensator C 110. In at least some examples, the parameters x of the predistorter are updated using an iterative updating procedure, for instance, using a stochastic gradient procedure. In other examples, the parameters x are updated in a batch procedure, for example, using a Gramian matrix inversion procedure.

In the iterative procedure, an objective is to optimize (i e, minimize) an objective

${E(x)} = {{\frac{\rho}{N}{x}^{2}} + {\frac{1}{N}{\sum\limits_{t \in T}{{{e\lbrack t\rbrack} - {\sum\limits_{k = 1}^{n}{x_{k}{B_{k}\left( {{q_{y}\lbrack t\rbrack},{p\lbrack t\rbrack}} \right)}}}}}^{2}}}}$ where e[t]=v[t]−y[t]

The sum is either over past time samples, or in the case of use of the buffer V, over the time samples represented in the buffer.

In a stochastic gradient approach, uniformly randomly selected times τ are selected and the current parameters values x are updated as x→x−α∇V _(τ)(x) where ∇_(τ)(x)=a[τ](a(τ)x−e[τ]) and a(τ)=[B ₁(q _(y)[t],p[t]), . . . ,B _(N)(q _(y)[t],p[t])]

The step size α may be adjusted (e.g., reduced) over time. Further details regarding the various optimization procedures that may be used in conjunction with the system 100 and the implementations of one or more of the agile adaptation techniques described herein are provided below.

Turning back to FIG. 1, the system 100 illustrated in FIG. 1 may be at least part of a digital front end of a device or system (such as a network node, e.g., a WWAN base station or a WLAN access point, or of a mobile device). As such, the system 100 may be electrically coupled to communication modules implementing wired communications with remote devices (e.g., via network interfaces for wired network connection such as through ethernet), or wireless communications with such remote devices. In some embodiments, the system 100 may be used locally within a device or system to perform processing (predistortion processing or otherwise) on various signals produced locally in order to generate a resultant processed signal (whether or not that resultant signal is then communicated to a remote device).

In some embodiments, at least some functionality of the linearization system 100 may be implemented using a controller (e.g., a processor-based controller) included with one or more of the modules of the system 100 (the compensator C 110, the adaptation section A 140, or any of the other modules of the system 100). The controller may be operatively coupled to the various modules or units of the system 100, and be configured to compute approximations of the gradient of the Error function (as will be discusses in greater detail below), update DPD coefficients based on those approximations, and perform digital predistortion using the continually updated DPD coefficients. The controller may include one or more microprocessors, microcontrollers, and/or digital signal processors that provide processing functionality, as well as other computation and control functionality. The controller may also include special purpose logic circuitry, e.g., an FPGA (field programmable gate array), an ASIC (application-specific integrated circuit), a DSP processor, a graphics processing unit (GPU), an accelerated processing unit (APU), an application processor, customized dedicated circuitry, etc., to implement, at least in part, the processes and functionality for the system 100. The controller may also include memory for storing data and software instructions for executing programmed functionality within the device. Generally speaking, a computer accessible storage medium may include any non-transitory storage media accessible by a computer during use to provide instructions and/or data to the computer. For example, a computer accessible storage medium may include storage media such as magnetic or optical disks and semiconductor (solid-state) memories, DRAM, SRAM, etc.

The approaches outlined can be used separately or in conjunction with any of the other DPD adaptation approaches discussed herein. Thus, for example, an embodiment of the approaches and techniques described herein include a method for digital predistortion comprising collecting data for a wireless device comprising a transmit chain that includes at least one power amplifier that produces output with non-linear distortions, the data including data records with each data record including input signal data and corresponding output signal data for the transmit chain and further being associated with a set of operating conditions. The each data record is associated with a variety metric representative of an extent of distinguishability of the respective data record from other collected data records. Collecting the data includes collecting the data in one or more spectral channels, the one or more spectral channels located within a wide band spectral range, with each of the one or more spectral channels having a respective bandwidth smaller than a range bandwidth of the spectral range. The method further includes selecting at least some of the data records according to a variety criterion computed based on the variety metric for the each data record, and computing, based on the selected at least some of the data records, digital predistortion parameters used to control a digital predistorter of the wireless device. Computing the digital predistortion parameters includes obtaining a first set of digital predistortion (DPD) basis functions for controlling operations of the digital predistorter, determining an adapted set of DPD basis functions, derived from the first set of DPD basis functions, that depend on operating conditions values corresponding to the selected at least some of the data records (for example, averages of corresponding operating condition values from each of the selected records may be computed), and configuring the digital predistorter with the digital predistortion parameters determined for the adapted set of DPD basis functions based, at least in part, on the selected at least some of the data records. The method additionally includes applying the digital predistorter, configured with the computed digital predistortion parameters, to subsequent input signals.

Additional details regarding the various agile adaptation techniques described herein follow.

First Technique: Wideband Adaptation Based on Narrow Band Measurements

As discussed herein, agile estimation of digital predistortion (DPD) parameters may be achieved with a universal parametrized digital predistorter, operable over a wideband spectral range, W, that can be implemented based on samples collected over portions of the wideband spectral range (e.g., based on samples collected in narrow subbands of the spectral range) at different time intervals. Such samples can be collected using an adjustable narrow band wireless communication device (configured to operate in different subbands within the spectral range).

Thus, in some embodiments, a method for digital predistortion of signals, provided to a radio frequency (RF) transmission path (e.g., the transmission section F 120) configured to transmit radio signals in a plurality of subbands within a spectral range, is disclosed. The method includes configuring a digital predistorter (such as the compensator C 110 of FIG. 1) for predistorting signals comprising arbitrary spectral content within the spectral range. The configuring includes acquiring one or more data samples representing operation of the RF transmission path to transmit radio signals in different subbands, with each data sample including a digital input signal provided to the RF transmission path for transmission, and a digital sensed (e.g., sensed by the sensing section S 130) signal representing a sensing of a radio transmission in response to the provided digital signal. Each digital input signal represents spectral content concentrated in a respective single corresponding subband within the spectral range (e.g., being concentrated or confined to a defined frequency range for a single communication channel or sub-carrier). Configuring the digital predistorter further includes updating parameters of the digital predistorter (e.g., through an adaptation procedure generally realized using the adaptation section A 140) according to the acquired one or more data samples to mitigate non-linear characteristics of the RF transmission path. The procedure further includes receiving a further input signal for transmission by the RF transmission path, the further input signal representing spectral content in a particular subband within the spectral range, and using the configured predistorter to process the further input signal to yield a predistorted signal for providing to the RF transmission path.

Thus, in implementations of the first technique, a universal parametrized predistorter is realized that is robust to changing spectral conditions. That is, in the approach of the first technique, a single predistorter implementation may be realized that can operate at different sub-carriers, and does not necessarily have to adapt to new spectral conditions as a channel's spectral characteristics (carrier frequency, bandwidth, communication protocols uses, average output power, etc.) change (for example, as a result of moving to a different cell, changing to a different band within the same cellular area due to changes caused by base station assignments, etc.) Because spectral characteristics can change frequently, using a robust universal digital predistorter configured to operate over the entire spectral range, W, can reduce computation complexity that would result from having to adapt to every change in spectral operating conditions.

To realize a universal parametrized digital predistorter, samples from subbands of the spectral range (i.e., bands within the spectral range, W, that have bandwidths narrower than the bandwidth W associated with the spectral range) are collected. The collection of samples can be performed at design time through a methodical traversal of the spectral range, or may be collected during normal operation of the wireless device (comprising the digital predistorter) as the wireless device hops from one subband to another, depending on the specific use of the wireless device, availability of spectral resources (e.g., whether some available bands are assigned to other wireless devices), and existing environmental conditions. When the wireless device (comprising a digital predistortion system such as the system 100 of FIG. 1) dynamically collects samples during normal operation, and occasionally (at regular or irregular intervals) updates DPD parameters, collected samples may be stored in a variety matrix such as the one described herein (i.e., to collect samples that promote sample diversity) along with operating conditions (that may include spectral operating conditions). Occasionally, the digital predistortion system uses the currently stored records to update the DPD coefficients. The updated coefficients would, in such examples, reflect a good variety of samples, including possible a variety of samples collected from diverse subbands. Alternatively, the wireless device may collect samples according to some other scheme that obtains representative samples obtained for multiple subbands.

With reference to FIG. 2, a diagram 200 illustrating collection of samples from narrow subbands (210, 212, and 214) of a wideband spectral range having a bandwidth of 100 MHz is provided. In this diagram, the wireless device hops (or systematically slides) into different subbands, and collects samples from the narrowband subbands. In addition to input samples (and/or observed samples), the operating conditions, which may include spectral information, can be collected and selectively stored (e.g., based on computed variety metrics, according to which a decision is made on whether particular samples should be stored or discarded). Occasionally, at regular or irregular intervals, the adaptation process uses the narrowband samples to generate or update the predistortion parameters, e.g., DPD coefficients that weigh the basis functions (which may be adaptive basis functions, as described herein). Sample data collection can be performed, in some embodiments, by sampling input (and/or observed output of the transmission chain) at full wideband bandwidth (i.e., at the Nyquist frequency for the spectral range W) or at the narrower bandwidth of the subband (optionally upsampling the collected samples during, for example, the DPD parameter updating procedure).

FIG. 3 includes graphs 310 and 320 comparing the performance results of DPD compensation systems that were implemented using two different approaches. In the first implementation, corresponding to the graph 310, the DPD coefficients were generated using wideband signal samples (i.e., the DPD system was trained using samples of a wideband signal). In the second implementation, corresponding to the graph 320, DPD coefficients were generated using samples collected from multiple narrow-band subbands. In the implementations corresponding to the graphs 310 and 320, derivation of the DPD coefficients was performed by accumulating and averaging new and old correlation matrices. As the two graphs 310 and 320 show, the DPD system trained using samples from multiple narrow-band subbands achieved a generally higher and more robust performance (based on the Adjacent Channel Leakage Ratio, or ACLR, behavior for the two implementations).

With reference next to FIG. 4, a flowchart of an example procedure 400 for digital predistortion of signals provided to a radio frequency (RF) transmission path, configured to transmit radio signals in a plurality of subbands within a spectral range, is shown. The procedure 400 includes configuring 410 a digital predistorter (such as the one realized by the system 100 of FIG. 1) for predistorting signals comprising arbitrary spectral content within the spectral range. The configuring operations 410 of the procedure 400 include acquiring 412 one or more data samples representing operation of the RF transmission path to transmit radio signals in different subbands. Each data sample includes a digital input signal (e.g., u[t] and/or v[t] shown in FIG. 1, and also illustrated in FIG. 4) provided to the RF transmission path for transmission and a digital sensed signal (e.g., y[t] shown in FIGS. 1 and 4) representing a sensing of a radio transmission in response to the provided digital signal. Each digital input signal represents spectral content concentrated in a respective single corresponding subband within the spectral range. The spectral information corresponding to the subband is shown in FIG. 4 as information that can optionally be provided for the purposes of the adaptation process. That spectral information can also be provided as part of the p[t] data depicted in FIG. 1. Thus, in some embodiments, acquiring the one or more data samples representing the operation of the RF transmission path in the different subbands may further include collecting for the one or more data samples corresponding operating conditions parameters, with the corresponding operating conditions parameters including one or more of, for example, a corresponding temperature of the RF transmission path, a corresponding average power, a corresponding power supply voltage (V_(DD)) associated with at least one power amplifier of the RF transmission path, a corresponding Voltage Standing Wave Ratio (VSWR) associated with the RF transmission path, age of the collected each of the sets of samples, and/or corresponding channel information in which the one or more data samples is collected. In some examples, the spectral range, W, includes n RF bands of approximately equal-sized bandwidths.

In some embodiments, acquiring the one or more data samples may include controllably adjusting a wireless device comprising the RF transmission path to operate in one or more of the different subbands within the spectral range (e.g., by controlling he modulation/demodulation functionality through adjustment of a local oscillator, adjusting bandpass filtering operations, etc.), and, at each of the one or more of the different subbands, measuring at least one set of samples. Acquiring the one or more data samples may include measuring during runtime, subsequent to an initial computation of the parameters of the digital predistorter, additional sets of samples at additional subbands within the spectral range, and updating the parameters of the digital predistorter further based on at least some of the additional sets of samples. In some embodiments, acquiring the one or more data samples representing operation of the RF transmission path to transmit radio signals in different subbands may include selecting the one or more data samples from a plurality of stored data samples according to a variety criterion computed based on the variety metric for the each data record, with the plurality of records being arranged in a variety database structure in which each record corresponds to one of the one or more data samples.

The configuring operation 410 also include updating 414 parameters of the digital predistorter (e.g., DPD coefficients that weigh basis functions) according to the acquired one or more data samples to mitigate non-linear characteristics of the RF transmission path. For example, in some embodiments, the basis functions are adaptive basis functions (basis functions that can be expanded from a base set, according to, for example, operating conditions) that include complex terms that depend on a τ-dimensional complex vector q corresponding to complex input samples to the digital distorter, and real terms that depend on a d-dimensional real vector p representative of values of operating conditions associated with operation of the RF transmission path. The value d is representative of a number of parameters of the operating conditions associated with the operation of the RF transmission path. In such embodiments, a total number of the basis functions used to digitally predistort input signals is proportional to the product of a number of terms, b, in a first set (i.e., a base set) of basis functions, and d+1, such that the total number of basis functions is proportional to b·(d+1).

With continued reference to FIG. 4, the procedure 400 further includes receiving 420 a further input signal (shown as the signal u_(subsequent)[t] in FIG. 4) for transmission by the RF transmission path, with the further input signal representing spectral content in a particular subband within the spectral range. The procedure 400 additionally includes using 430 the configured predistorter (configured according to the operations 410) to process the further input signal to yield a predistorted signal (shown in FIG. 4 as v_(subsequent)[t]) for providing to the RF transmission path.

In some examples, using the configured predistorter to process the further input signal may include predistorting the further input signal according to:

${{v\lbrack t\rbrack} = {{u\lbrack t\rbrack} + {\sum\limits_{k = 1}^{n}{x_{k}{B_{k}\left( {{q_{u}\lbrack t\rbrack},{p\lbrack t\rbrack}} \right)}}}}},{{{where}\mspace{14mu}{q_{u}\lbrack t\rbrack}} = \begin{bmatrix} {u\left\lbrack t + l - 1 \right.} \\ \vdots \\ {u\left\lbrack t + l - \tau \right.} \end{bmatrix}}$ where v[t] is an output of the digital predistorter at time t, u[t] is an input sample to the digital predistorter at the time t, B_(k) are the basis functions for k=1 to n, q_(u)[t] are stacks of input samples, p[t] are d-dimensional real vectors of operating conditions (e.g., environmental characteristics) at the time t, and x_(k) are the computed digital predistortion parameters to weigh the basis functions.

Updating the parameters of the digital predisorter may include deriving, according to an optimization process that is based, at least in part, on the acquired one or more data samples measured in the different subbands, the digital predistortion parameters, x_(k), to weigh the basis functions B_(k). In some examples, deriving according to the optimization process the parameters, x_(k), of the digital predistorter may include adaptively adjusting the digital predistortion parameters, x_(k), through minimization of a regularized mean square error function E(x) according to:

${{E(x)} = {{\frac{\rho}{n}{x}^{2}} + {\frac{1}{N}{\sum\limits_{t \in T}{{{v\lbrack t\rbrack} - {y\lbrack t\rbrack} - {\sum\limits_{k = 1}^{n}{x_{k}{B_{k}\left( {{q_{u}\lbrack t\rbrack},{p\lbrack t\rbrack}} \right)}}}}}^{2}}}}},$ where ρ>0 is a regularization factor, v[t] is an output of the digital predistorter at the time t, and y[t] is the output of an observation receiver coupled to an output of the RF transmission path.

Adaptively adjusting the digital predistortion parameters, x_(k), may include, for example, deriving correlation matrices G and L, when n<<100, with n being representative of the number of basis functions, according to:

${G = {\frac{1}{N}{\sum\limits_{t \in T}{{B\lbrack t\rbrack}{B\lbrack t\rbrack}^{\prime}}}}},{L = {\frac{1}{N}{\sum\limits_{t \in T}{\left( {{v\lbrack t\rbrack} - {y\lbrack t\rbrack}} \right)^{\prime}{B\lbrack t\rbrack}}}}},{{{where}\mspace{14mu}{B\lbrack t\rbrack}} = \begin{bmatrix} {B_{1}\left( {{q_{y}\lbrack t\rbrack},{p\lbrack t\rbrack}} \right)} \\ \vdots \\ {B_{n}\left( {{q_{y}\lbrack t\rbrack},{p\lbrack t\rbrack}} \right)} \end{bmatrix}},$ and deriving the digital predistortion parameters, x_(k), according to

$\overset{\_}{x} = {\left( {{\left( \frac{p}{n} \right)I} + G} \right)^{- 1}{L.}}$ In the above expression, x is a column vector of the complex conjugated x_(k).

Adaptively adjusting the digital predistortion parameters, x_(k), may include, in some embodiments, approximating the digital predistortion parameters, x_(k), according to a stochastic gradient process.

Second Technique: Parameterized Digital Predistortion Basis Functions

As noted, the coefficients derived by the adaptation section A 140 are used to weigh the basis functions used by compensator (predistorter) C 110 to predistort the input u. Predistorting the signal u, by the compensator C 110 requires, in some embodiments, minimizing the difference between y (the output of the sensing section S 130) and u (the input samples) as small as possible. To add agility and versatility to the predistortion operations, some example embodiments described herein utilize a set of adaptable basis functions to build an expanded set of functions based on contemporaneous operating conditions. In other words, as the operating conditions (e.g., environmental characteristics) change, the particular set of basis functions operating on the input samples is adjusted in accordance with the existing operating conditions.

An example approach that uses an adapted set of basis functions that depends on operating conditions (to achieve agile estimation of DPD parameters, improve responsiveness speed to changing conditions, and increase robustness of the DPD system) is based on deriving an adapted set of basis functions from an initial/first set of DPD basis functions. The first set of DPD functions may be a pre-determined, carefully selected set of functions B_(k)(•) that is dependent on the input samples to the DPD system (and is configured to provide good DPD performance, e.g., for a particular transmit chain or wireless device), but is independent of operating conditions. Thus, the initial (base) set of basis functions (from which descendent/derivative sets of basis functions can be spawned) comprises terms that are applied to input samples to the DPD system (i.e., u[t−k], . . . , u[t−1], u[t], u[t+1], . . . ). Spin-off (child or derivative) basis functions can then be constructed using the initial set of functions, with added terms that are functions of the operating conditions (e.g., parameters representative of environmental characteristics). As a result, the size of adaptive sets of basis functions increases as more terms dependent on operating conditions are added to the initial (base) set of basis functions.

A common representation of the compensator C 110 (of FIG. 1) is given by:

${v\lbrack t\rbrack} = {{u\lbrack t\rbrack} + {\sum\limits_{k = 1}^{n}{x_{k}{B_{k}\left( {{q_{u}\lbrack t\rbrack},{p\lbrack t\rbrack}} \right)}}}}$ where

${q_{u}\lbrack t\rbrack} = \begin{bmatrix} {u\lbrack t + l - 1\rbrack} \\ {u\lbrack t + l - 2\rbrack} \\ \vdots \\ {u\lbrack t + l - \tau\rbrack} \end{bmatrix}$ is the stack of recent (around t) baseband input samples. As noted above, the vector p[t] is d-dimensional real vector of operating conditions at time t, x_(k) (k=1, . . . , n) are the elements of a complex n-dimensional vector x∈C^(n) (the coefficients of the compensator), and B_(k)(q,p) are carefully selected “basis functions” of τ-dimensional complex vector q and d-dimensional real vector p (a “careful” selection means that the operation of computing a single sample v[t] of the corrected signal takes significantly less than n scalar multiplications; examples of carefully selected basis functions are provided in PCT Application No. PCT/US2019/031714, filed on May 10, 2019, titled “Digital Compensation for a Non-Linear System”). Two example operating conditions (i.e., in a situations where d=2) that may be used to form the expanded set of basis functions may include the drain voltage (V_(DD)) associated with the at least one power amplifier of the wireless device, and the temperature, in which case the vector p[t] may be represented as follows:

${p\lbrack t\rbrack} = {\begin{bmatrix} {p_{1}\lbrack t\rbrack} \\ {p_{2}\lbrack t\rbrack} \end{bmatrix} = \begin{bmatrix} {V_{DD}\lbrack t\rbrack} \\ {Tem{p\lbrack t\rbrack}} \end{bmatrix}}$

In some examples, the formation of the adapted, spinoff, set of basis functions (comprising b·(d+1) basis functions, where b is the total number of terms in the initial/base set of basis functions, and d is the number of operating conditions that are used to generate the expanded, spinoff, set of basis functions) is achieved by forming one group of basis functions (comprising b terms) as the unmodified initial (“base”) set of basis functions, and forming d additional groups of basis functions (each with b basis function terms) as the product of the respective measured operating conditions and the initial/base set of basis functions. Thus, an expanded set of adaptive basis functions, B_(k)(q_(u)[t], p[t]), for a situation involving use of two (2) operating conditions (e.g., p₁ and p₂, corresponding to V_(DD) and temperature) can be represented as: B _(3(k−1)+1)(q)=B _(k0)(q) B _(3(k−1)+2)(q)=B _(k0)(q)·p ₁ B _(3(k−1)+3)(q)=B _(k0)(q)·p ₂

Other ways to generate an expanded set of basis functions, e.g., using other formulations and transformations to express as functions (linear or non-linear) of the operating conditions (p₁, p₂, . . . , p_(n)), with such formulations and transformations being applied to the initial/base set of basis functions, may be used instead of or in addition to the above formulation involving the product of the operating conditions of p₁ and p₂ with the base set of basis functions. Furthermore, some groups of expanded basis function may be formed as a function of multiple operating conditions (e.g., each group of the expanded set may be expanded as a function/transformation involving combinations of p₁, p₂, p_(n)). It is to be noted that expansion of basis functions works best when there is some linear dependency between the operating conditions (e.g., p₁ and p₂, in the above example).

In operation, a compensator, such as the compensator C 110, operating according to an inverse model of the nonlinear distortion of the transmit chain (e.g., realized by the block F 120 of FIG. 1), and using adaptive basis functions (expanded and adapted according to contemporaneous operating condition values) that are weighed by regularly updated DPD coefficients, receives the input signal, u, and generates the intermediate input signal, v as follows:

$v = {{2u} + {\sum\limits_{i = 1}^{n}{x_{i}{f_{i}(u)}}}}$ where f_(i)(•) is the i^(th) basis function of the basis functions, and x_(i) is the i^(th) parameter (e.g., the i^(th) weight) corresponding to the i^(th) basis function. Each of the basis functions (at least in the base set) is a linear function (e.g., u(t−1)) or a non-linear function (e.g., |u(t)|²) of the input, u, which may include memory (e.g., u(t)*u(t−1)).

In some implementations, the adaptation process implemented by the system 100 of FIG. 1 adjusts the coefficients set of coefficients, x_(k) (where k represents a time instance or interval during which that set x_(k) is being used) of the predistorter by observing the measurement samples y[t] and trying to minimize the regularized mean square error according to:

${{E(x)} = {{\frac{\rho}{N}{x}^{2}} + {\frac{1}{N}{\sum\limits_{t \in T}{{{e\lbrack t\rbrack} - {\sum\limits_{k = 1}^{n}{x_{k}{B_{k}\left( {{q_{y}\lbrack t\rbrack},{p\lbrack t\rbrack}} \right)}}}}}^{2}}}}},$ where

${{q_{y}\lbrack t\rbrack} = \begin{bmatrix} {y\lbrack t + l - 1\rbrack} \\ {y\lbrack t + l - 2\rbrack} \\ \vdots \\ {y\lbrack t + l - \tau\rbrack} \end{bmatrix}},$

In the above expression, e[t]=v[t]−y[t] is the input/output mismatch in the observed data, ρ>0 is the regularization factor (used to prevent poor conditioning of the optimization problem, which tends to result in unstable interaction between the adaptation process and the compensator application, as well as unreasonably large entries of the optimal x+=argmin E(x)), T={t} is a set of time samples, and N is the number of elements in T.

When n is small (e.g., n<<100), minimization of E(x) can be carried out by computing the following correlation matrices:

${G = {\frac{1}{N}{\sum\limits_{t \in T}{{B\lbrack t\rbrack}{B\lbrack t\rbrack}^{\prime}}}}},{L = {\frac{1}{N}{\sum\limits_{t \in T}{\left( {{v\lbrack t\rbrack} - {y\lbrack t\rbrack}} \right)^{\prime}{B\lbrack t\rbrack}}}}},$ where B[t] is a vector of the basis functions (i.e., a vector of complex entries, with one entry per basis function of n basis functions:

${B\lbrack t\rbrack} = \begin{bmatrix} {B_{1}\left( q_{y}\lbrack t\rbrack,p\lbrack t\rbrack \right)} \\ \vdots \\ {B_{n}\left( {{q_{y}\lbrack t\rbrack},{p\lbrack t\rbrack}} \right)} \end{bmatrix}$

In some example, the DPD coefficients to weigh the basis functions can then be computed according to the following expression, which includes a term dependent on ρ to regularize the optimization:

$\overset{\sim}{x} = {\left( {{\left( \frac{p}{n} \right)I} + G} \right)^{- 1}L}$

Alternatively, in some embodiments (e.g., when n is large), the DPD coefficients to weigh the adaptive basis functions can be computed according to a stochastic gradient process. An example stochastic gradient process that can be used is as follows.

The function E(x) to be minimized can be expressed in the form:

${{E(x)} = {{{\rho{x}^{2}} + {\frac{1}{N}{\sum\limits_{t \in T}{{{e\lbrack t\rbrack} - {{a\lbrack t\rbrack}{x\lbrack t\rbrack}}}}^{2}}}} = {R + {2{{Re}\left\lbrack {L^{\prime}x} \right\rbrack}} + {x^{\prime}{Gx}} + {\rho{x}^{2}}}}},$ where a[t]=[B₁(q_(y)[t],p[t]), B₂(q_(y)[t],p[t]), . . . , B_(n)(q_(y)[t],p[t])] are the complex 1-by-n matrices formed by the basis function samples, and

${G = {\frac{1}{N}{\sum{{a\lbrack t\rbrack}^{\prime}{a\lbrack t\rbrack}}}}},{L = {\frac{1}{N}{\sum{{a\lbrack t\rbrack}^{\prime}{e\lbrack t\rbrack}}}}},{R = {\frac{1}{N}{\sum{{e\lbrack t\rbrack}}^{2}}}}$ are the components of the Gramian matrix computed from the data.

As time progresses, the data represented by the samples (q_(y)[t],p[t],e[t]), t ∈T, gets refreshed with the outcomes of new measurements. In some situations, this can be done by either re-computing the Gramian, by averaging new and old Gramians, or by replacing the old (q_(y)[t],p[t],e[t]) samples with the new ones.

In some example implementations, the adaptation process used for a DPD implementation attempts to find the vector x of compensator coefficients as an approximation of the optimal x_(*)=argmin E(x). An approach to finding x_(*)=argmin E(x) is by computing it as the (unique) solution x_(*)=(pI_(n)+G)⁻¹L of the linear equation (pI_(n)+G)x=L, where I_(n) denotes the n-by-n identity matrix. This approach requires computation (and storage) of the n-by-n Hermitian matrix G, which makes it impractical even for very simple compensators (for example, a compensator with n=100 complex coefficients will require computing and storing n²=10000 real coefficients of G). In addition, solving large systems of linear equations on FPGA or ASIC is a challenge on its own, and does not resolve the challenge of large storage required for Gramians.

Alternatively, the gradient of E=E(x) at a point x∈C^(n) is given by:

${{\nabla{E(x)}} = {{2\left( {{\rho x} + {Gx} - L} \right)} = {2\left\lbrack {{\rho x} + {\frac{1}{N}{\sum\limits_{t \in T}{{a\lbrack t\rbrack}^{\prime}\left( {{{a\lbrack t\rbrack}x} - {e\lbrack t\rbrack}} \right)}}}} \right\rbrack}}},$ and can be computed and stored as an n-by-1 complex matrix, which reduces the storage requirements, compared to those of the Gramian.

Once computed, the gradient can be used to adjust the current guess {tilde over (x)} at the optimal x, according to {tilde over (x)}⁺={tilde over (x)}−γ∇E({tilde over (x)}), making it slightly better (as long as the step size parameter γ is sufficiently small). With the updated guess {tilde over (x)}={tilde over (x)}⁺, the gradient can be re-computed, and another adjustment can be made. Repeating the process sufficiently many times leads to an arbitrarily accurate approximation {tilde over (x)} of x_(*). Thus, in approaches based on a gradient of the error function E(x), a new/subsequent set of coefficients to weigh the basis functions applied by the compensator (e.g., the compensator C 110) can be estimated/approximated as a change of the current coefficients by a pre-determined multiple (i.e., a step-size γ) of the value of the gradient of E(x).

Re-computing the gradient requires a significant number of operations each time. Since the number of steps of gradient search, before it reaches a satisfactory level of optimality, tends to be significant, compensator adaptation based on gradient search may be slow. Accordingly, in some embodiments, to reduce the amount of computation required to derive a current value of the gradient of E(x) (based on which updated coefficients are to be computed), a sub-gradient approach to minimize E(x) may be used to approximate the DPD coefficients x. Particularly, a stochastic gradient approach to minimize the function E=E(x) is implemented based on the fact that the gradient ∇E(x) is the expected value of the vector random variable ξ=2[ρx+a[τ]′(a[τ]x−e[τ])], where τ is the random variable uniformly distributed on the finite set T={t}, and that random samples of ξ (which are much cheaper, from a computational effort, to compute than the true gradient ∇E(x)) can replace ∇E(x) in the gradient search iterations without a significant loss of convergence rate. In the above expression, a[τ] represents a matrix of values resulting from application of the basis functions, as weighed by the DPD coefficients x, to the input sample(s) at a randomly selected instant corresponding to the random variable τ, and e[τ] is the input/output mismatch in the observed data between an output of an observer receiver and the output of the digital compensator (e.g., the output y of the sensing section S 130, and the output v of the compensator C 110).

Thus, a stochastic gradient approach that may be used includes constructing a sequence of random values {tilde over (x)}₀, {tilde over (x)}₁ . . . (taking values in C^(n)), and updating the current coefficients at an instance k (namely, {tilde over (x)}_(k)) as follows: {tilde over (x)} _(k+1) ={tilde over (x)} _(k) +αa[τ_(k)]′(e[τ_(k)]−a[τ_(k)]{tilde over (x)} _(k))−αρ{tilde over (x)} _(k) ,{tilde over (x)} ₀=0, where τ₁, τ₂, . . . , τ_(n) are independent random variables uniformly distributed on T={t}, ρ>0 is the regularization constant from the definition of E=E(x), and α>0 is a multiple (which may be a constant value that can be adjusted at regular or irregular periods) such that α(ρ+|a[t]|²)<2 for every t ∈T. The value α represents the step size that is applied to a value based on the expected value of the vector random variable used to approximate the gradient ∇E(x).

The expected value x _(k)=E[{tilde over (x)}_(k)] can be shown to converge to x_(*)=argmin E(x) as k→∞. An averaging operation may be used according to the expression: {tilde over (y)} _(k+1) ={tilde over (y)} _(k)+ϵ({tilde over (x)} _(k) −{tilde over (y)} _(k)), with ϵ∈(0,1].

The difference between {tilde over (y)}_(k) and x_(*) is small for large k, as long as ϵ>0 is small enough. This approach to minimizing E(x) is referred to as a “projection” method, since the map x→x+|a[t]|⁻² a[t]′(e[t]−a[t]x[t]) projects x onto the hyperplane a[t]x=e[t].

In practical implementations of the processes to compute approximation of the gradient of E(x), and to thus update the current coefficients, τ_(k) are generated as a pseudo-random sequence of samples, and the calculations of {tilde over (y)}_(k) can be eliminated (which formally corresponds to ϵ=1, i.e., {tilde over (y)}_(k)={tilde over (x)}_(k−1)). As a rule, this requires using a value of α with smaller minimal upper bound for α(ρ+|a[t]|²) (for example, α(ρ+|a[t]|²)<1, or α(ρ+|a[t]|²)<0.5). In some embodiments, the values of α is sometime adjusted, depending on the progress made by the stochastic gradient optimization process, where the progress is measured by comparing the average values of |e[τ_(k)]| and |e[τ_(k)]−a[τ_(k)]{tilde over (x)}_(k)|. In some embodiments, a regular update of the set of the optimization problem parameters a[t], e[t] may be performed to replace the data samples a[t], e[t] observed in the past with new observations.

Because, in the stochastic gradient approach described herein, the DPD coefficients are dynamically updated (reflecting the fact that the behavior of the system 100 changes over time), the expected values for the vector random value ξ represents a good approximation or estimate for the gradient ∇E(x), thus allowing for cheaper computation (from a computational effort perspective) of the values used to update the DPD coefficient x used by the compensator of the system 100.

In some embodiments, the value of the step size (that multiplies the approximation determined for the gradient ∇E(x)) may be adjusted based on a measure of progress of the stochastic gradient process, with the measure of progress being determined by, for example, comparing an average values of |e[τ_(k)]| and |e[τ_(k)]−a[τ_(k)]{tilde over (x)}_(k)|. In some examples, the process starts off with a relatively larger step size and then gradually reduces the step size.

With reference next to FIG. 5, a flowchart of an example procedure 500 for digital predistortion, based on the use of adaptive basis functions (e.g., formed using operating conditions of a digital predistortion system performing the procedure 500) is shown. The procedure 500 includes obtaining 510 a first set of digital predistortion (DPD) basis functions (shown as the base basis function B_(k)(q), which is expressed as a function of the samples of the input signal u, namely, q_(u)[t]=(u[t+l−1], u[t+l−2], . . . , u[t+l−τ])) for controlling operations of a digital predistorter, of a wireless device, operating on a received at least one input signal. The at least one input signal is directed to a power amplification system (such as the system 100 of FIG. 1) comprising a transmit chain (such as the transmission section F 120 of FIG. 1) with at least one power amplifier that produces output with non-linear distortions.

With continued reference to FIG. 5, the procedure 500 additionally includes determining 520 an adapted set of DPD basis functions, derived from the first set of DPD basis functions, that depends on operating conditions (shown as being provided through the data p[t] received at the block 520) associated with operation of the wireless device, for controlling operation of the digital predistorter. In some examples, the operating conditions may include one or more of, for example, channel information associated with a corresponding channel used by the wireless device to process the at least one input signal, a corresponding temperature, a corresponding average power, a corresponding drain voltage (V_(DD)) associated with the at least one power amplifier of the wireless device, a corresponding Voltage Standing Wave Ratio (VSWR) associated with the wireless device, and/or age of collected samples used for determining the DPD coefficients.

As discussed herein, the adaptive set of basis function if formed through transforming an original set of carefully selected basis functions that includes complex terms that depend on a τ-dimensional complex vector q (corresponding to complex input samples to the digital distorter) to an expanded set of functions through transformations that are implemented, for example, by real terms that depend on a d-dimensional real vector p representative of values of the operating conditions (including environmental characteristics) associated with the operation of the wireless device on the received at least one input signal. The value d is representative of a number of parameters of the operating conditions, associated with the operation of the wireless device, that are being used to expand the base set of the basis functions. In such embodiments, a total number of basis functions (in the expanded set) used to digitally predistort the at least one input signal is proportional to the product of a number of terms, b, in the first set of basis functions, and d+1, such that the total number of basis functions is proportional to b·(d+1). The basis functions thus comprise (d+1) groups of basis functions, with a first group of basis functions including, in some examples, the first set of basis functions without any modification, and with each of the remaining d groups of basis functions determined as the product of the first set of basis functions and one of the operating conditions (e.g., environmental characteristics) associated with the wireless device. In some other examples, other types of transformations, expressed as linear or non-linear functions of the operating conditions, may be used to expand the base set of basis functions.

With continued reference to FIG. 5, the procedure 500 also includes configuring the digital predistorter with DPD coefficients (illustrated as the output x_(k) of the block 530) determined for the adapted set of DPD basis functions based, at least in part, on observed samples (e.g., y[t] in FIG. 1, and as illustrated as an input to the block 530, and/or the input signal u[t] or the output v[t] of the compensator C 110 of FIG. 1) of the transmit chain responsive to the at least one input signal.

In some embodiments, configuring the predistorter with the DPD coefficients may include computing the DPD coefficients according to an optimization process that yields the DPD coefficients, x_(k), as weights to the basis functions that minimizes an error function based on the observed samples of the transmit chain and samples of the at least one input signal. In such embodiments, computing, according to the optimization process, the DPD coefficients, x_(k), may include adaptively adjusting the DPD coefficients, x_(k), through minimization of a regularized mean square error function E(x) according to:

${{E(x)} = {{\frac{\rho}{n}{x}^{2}} + {\frac{1}{N}{\sum\limits_{t \in T}{{{v\lbrack t\rbrack} - {y\lbrack t\rbrack} - {\sum\limits_{k = 1}^{n}{x_{k}{B_{k}\left( {{q_{y}\lbrack t\rbrack},{p\lbrack t\rbrack}} \right)}}}}}^{2}}}}},{{{where}{\mspace{14mu}\;}{q_{y}\lbrack t\rbrack}} = \begin{bmatrix} {y\lbrack t + l - 1\rbrack} \\ \vdots \\ {y\lbrack t + l - \tau\rbrack} \end{bmatrix}}$ where ρ>0 is a regularization factor, B_(k) are the basis functions for k=1 to n, q_(y) [t] are stacks of input samples at time t, p[t] are d-dimensional real vectors of the operating conditions (which may include environmental characteristics) at the time t, v[t] is an output of a digital predistortion (DPD) compensator at the time t, and y[t] is the output of an observation receiver coupled to an output of the transmit chain.

As described herein, adaptively adjusting the DPD coefficients, x_(k), may include deriving correlation matrices G and L, when n<<100, with n being representative of the number of basis functions, according to:

${G = {\frac{1}{N}{\sum\limits_{t \in T}{{B\lbrack t\rbrack}{B\lbrack t\rbrack}^{\prime}}}}},{L = {\frac{1}{N}{\sum\limits_{t \in T}{\left( {{v\lbrack t\rbrack} - {y\lbrack t\rbrack}} \right)^{\prime}{B\lbrack t\rbrack}}}}},{{{where}\mspace{14mu}{B\lbrack t\rbrack}} = \begin{bmatrix} {B_{1}\left( {{q_{y}\lbrack t\rbrack},{p\lbrack t\rbrack}} \right)} \\ \vdots \\ {B_{1}\left( q_{y}\lbrack t\rbrack,p\lbrack t\rbrack \right)} \end{bmatrix}},$ and deriving the DPD coefficients, x_(k), according to

$\overset{\_}{x} = {\left( {{\left( \frac{p}{n} \right)I} + G} \right)^{- 1}{L.}}$

Alternatively, adaptively adjusting the DPD coefficients, x_(k), may include approximating the DPD coefficients, x_(k), according to a stochastic gradient process.

Third Technique: Variety-Based Adaptation Approach

As discussed herein, another technique/approach to achieve agile estimation of predistortion parameters, based on sparse data acquisition/availability, and/or to counter rapidly changing operating conditions, is through the implementation of a variety matrix that is used to maintain a limited number of sample points based on which adaptation processes are performed. To facilitate robust digital predistortion performance, the variety matrix techniques described herein use optimization data that includes input/output samples collected from a diverse set of operating conditions. The diversity (i.e., variety) of the data used to perform the optimization process (e.g., using one or more of the optimization processes described herein) is achieved using a diversity matrix that allows the collection of data samples that provide a greater diversity of the data set, and preserves historical data samples that have high relevance and usability for updating the DPD coefficient (samples' relevance can decrease over time, as the previously collected samples age).

Briefly, under the variety matrix approach, a DPD system, such as the system 100 depicted in FIG. 1, collects data (e.g., maintained in a memory storage device coupled to the adaptation section A 140, or to any other of the subsystems of FIG. 1) that includes data records, with each data record including input signal data and corresponding output signal data for the transmit chain (e.g., implemented via the transmission section F 120). The records may be maintained in a matrix-like data structure, and each may be associated with a set of operating conditions values (e.g., temperature, the power supply voltage associated with the power amplifier of the transmission section, load conditions, band attributes associated with the spectral band through which the data was collected or is to be transmitted to, a VSWR value, sampling rate, modulation mode, age of the samples, etc.), which may also be stored within the matrix. Each data record may be associated with a variety measure/score that is representative of an extent of distinguishability of the respective data record from other collected data records. Such a measure may be stored as a field or entry of the data records (e.g., within the matrix storing others portions of the collected data), or may be stored as a separate data structure, with the positions of different measures within the data structure corresponding to the particular records with which the different measures are respectively associated.

To manage the computational complexity involved with the optimization process used to compute adaptive DPD coefficients (that weigh the basis functions, such as the adaptive basis functions described herein), a subset of collected samples may need to be selected. This selection may be required in situations where the method of optimization is one based on use of correlation matrices (which, as was noted above, is computationally manageable when the number of data samples used in the optimization process is relatively small, e.g., n<<1000). Thus, in the variety-based approach discussed herein, the subset of records used for optimization is selected according to the variety measures associated with the data records. For example, in embodiments in which the variety measures are associated with the extent/degree of distinguishability of records from each other, the selection of the records to be used in the optimization can be one where the k records (from the available n records, where k<n) that provide the most distinguishability (i.e., the most diversity or variety) are selected. In some implementations, the selection of records is performed dynamically, concomitantly with receipt of data, by deciding whether to replace an existing record within a matrix structure, associated with a particular variety measure/score, with a newly arrived record whose variety measure indicates a higher degree of distinguishability than the existing record. In that situation, the previously stored record is replaced in favor of the new record. Further details about the procedure to replace records (to thus maintain records that increase or promote the extent of record diversity) are provided below.

Computation of the variety measure (which may be performed by the adaptation section A 140) may include computing, for the each data record, the variety measure based on distance metric representative of distances between data elements for a particular signal instance (e.g., the values representative of the samples and/or the values representative of operating conditions for the data record for the signal instance), and respective data elements of records (stored in the matrix) for other signal instances. The distance may be computed according to one or more formulation used for determining distances in a mathematical hyperspace of dimension n. Other processes for computing a variety measure/score may also be used. The computed variety measures can then be ranked, and, in implementations in which only a certain number (k) of the available signal instances are used for deriving the DPD coefficients, the k records with the highest variety measures/scores may be selected for use in the optimization process. Other criteria for selecting the k records, according to the variety measures associated with the data records may be used.

As noted, the collected data records, comprising input and output samples, and corresponding operating conditions, can be arranged in a matrix data structure, also referred to as a variety database V, as follows:

$V = \begin{bmatrix} q_{1} & q_{2} & \ldots & q_{m} \\ p_{1} & p_{2} & \ldots & p_{m} \\ e_{1} & e_{2} & \ldots & e_{m} \end{bmatrix}$ where q_(i)=q_(y)[t_(i)] is a complex vector of input samples collected for a signal i (in some examples, 5, 10, or any number of samples may be collected for every signal instance i), p_(i)=p[t_(i)] is the corresponding vector of operating conditions associated with the record corresponding to the signal instance i, and e_(i)=v[t_(i)]−y[t_(i)] is the difference between the resultant predistorted output of the compensator (such as the compensator C 110 of FIG. 1) and the observed output sample of the transmit chain (the e_(i) values are used for the optimization process to compute the adaptive DPD coefficients to weight the basis functions applied to input signals u[t] provided to the input of the compensator). The variety matrix V can be of any size, but in some typical examples the number of columns, m, is less than a 1000, which allows the optimization process for the adaptation procedure to be performed, for example, through the correlation matrices techniques described herein. As noted, if the number of distinct sample records in the database (distinct i) is relatively large, the computations complexity for optimizing the DPD coefficients using the correlation matrices technique may be too high, and, under such circumstances, a stochastic gradient process may be employed. The following V matrix example stored for each signal instance 3 samples, and two operating condition parameters (e.g., VDD and temperature), resulting in a matrix defining a 5-dimensional space.

$V = \begin{bmatrix} {u_{1}\left\lbrack {t + 1} \right\rbrack} & {u_{2}\left\lbrack {t + 1} \right\rbrack} & \ldots & {u_{m}\left\lbrack {t + 1} \right\rbrack} \\ {u_{1}\left\lbrack {t + 0} \right\rbrack} & {u_{2}\left\lbrack {t + 0} \right\rbrack} & \ldots & {u_{m}\left\lbrack {t + 0} \right\rbrack} \\ {u_{1}\left\lbrack {t - 1} \right\rbrack} & {u_{2}\left\lbrack {t - 1} \right\rbrack} & \ldots & {u_{m}\left\lbrack {t - 1} \right\rbrack} \\ {p_{1}(1)} & {p_{2}(1)} & \ldots & {p_{m}(1)} \\ {p_{1}(2)} & {p_{2}(2)} & \ldots & {p_{m}(2)} \\ \ldots & \ldots & \ldots & \ldots \end{bmatrix}$

Generally, however, each signal instance (stored in each of the columns of the V matrix) will include more than three (3) signal sample points in order to improve the accuracy and robustness of the optimization process to derive the DPD coefficients, x_(k).

In some embodiments, the matrix V may be accompanied by a minimal distance vector R, namely R=[r₁, r₂, . . . , r_(m)] which maintains variety measures (e.g., one measure for each of the columns/records of the variety matrix V) that can be computed as a distance function between various elements of a particular signal instance (e.g., a newly arrived signal instance), and any of the other columns/records of the variety matrix V. In such embodiments, the greater the distance between any two records, the greater the diversity score (i.e., the distinguishability) corresponding to the records used to optimize the DPD coefficients. In some embodiments, the minimal distance R vector entries may be combined with the V matrix (e.g., as an added row of the matrix V). Large variety measure values, as reflected by the vector R, may be indicative of a high degree of distinguishability between the corresponding columns (which correspond to the signal instances) and thus indicative of a good distribution of the samples over the possible space of operating conditions. Large variety can reduce sampling bias (or conversely increase the randomization of samples) that may have been caused if the sample instances collected are too closely lumped to each other (i.e., their corresponding operating condition are too closely lumped, which would result in a less robust predistorter).

As noted, in some embodiments, a matrix of a predetermined size (e.g., with k columns, and m rows, where each column represents a separate signal instance, and each row represents values associated with a signal instance, including the sample values and values representative of operating conditions) is initially populated with an initial data set (e.g., all ‘0’s), and gradually, as new samples arrive, their distance (according to one or more distance formulation) to data entries stored in the V matrix are computed, and a determination of whether the newly arrived sample offers more diversity (relative to currently stored sampling instances) is made. If the newly arrived sample would cause a higher diversity of data maintained in the V matrix, that newly arrived data replaces one of the columns in the matrix V.

An example of a process to determine (according to one or more criteria that are based on computed variety measures) if existing data within the matrix V is to be replaced by data corresponding to a newly arrived data entry (q_(y)[t], p[t],v[t]−y[t]) is as follows.

-   -   1) Compute h_(i)=dist_(*)((q_(i), p_(i)), (q_(y)[t],p[t])) for         i=1, . . . , m, where the approximate distance dist_(*)((q, p),         ({tilde over (q)}, {tilde over (p)})) denotes an easy-to-compute         measure of distance between the pairs (q, p) and ({tilde over         (q)}, {tilde over (p)}). In some embodiments, h_(i) may be         derived as a maximal absolute value of difference between the         corresponding individual components of q, p, {tilde over (q)},         {tilde over (p)} in a quantized grid.     -   2) Replace each entry r_(i) of R with min{r_(i), h_(i)}.     -   3) Find the index k ϵ{1, 2, . . . , m} of a minimal entry r_(k)         in the updated R. In this operation, the replacement process         identifies an entry from the updated R vector with the smallest         minimal distance. This entry will correspond to a column whose         relative distinguishability to other sample data points is low         (i.e., it is associated with a low variety measure).     -   4) Replace the k-th column of V matrix with a new column,         according to q_(k)=q_(y) [t], p_(k)=p[t], e_(k)=v[t]−y[t].     -   5) Replace the k-th entry of R with h_(i).

As noted, in some embodiments, one example of an operating condition parameter that may be used is the age of the existing data sample. In such embodiments, a relatively old existing sample data instance in the matrix V may result in a large distance value (i.e., a large h value) between that old existing sample data instance and the new data sample point. Thus, if the new data sample point replaces the old existing sample data instance in the matrix V, the corresponding entry in the R vector will have a relatively large h value, and consequently that just replaced column will be less likely to be replaced by the next sample data instance collected (i.e., because the corresponding entry is less likely to be the minimal entry in the R vector).

In some implementations, a decision may be made not to replace any of the existing columns with a new sampling data instance. For example, if it is determined that the computed variety measure for the new sampling data instance would result in a less diverse set of sample instances (e.g., because there is already a large group of sampling data points that are similar to the newly considered sampling data point), the replacement procedure may determine to retain all the existing records (e.g., matrix columns) and discard the newly collected sampling data instance.

Other approaches and techniques to control (e.g., replace) the content of the V matrix (by applying different variety criteria that are based on computed variety measures) may be implemented. Such approaches, including the above approach to replace columns of the V matrix, are configured to allow the most diverse data to be kept in the V matrix.

At regular or irregular time intervals, (e.g., every second, 10 seconds, 1 minute, 1 hour, etc.), an adaptation process such as the various processes discussed herein may be performed using the available sampling data points (including new sampling data points that have been added to the V matrix since the previous update). Alternatively, in some examples the adaptation process may be run after certain number of updates has been performed on the V matrix (e.g., after 50% of the columns of the V matrix have been replaced since the previously run adaptation process). As discussed herein, the adaptation process to determine DPD coefficients to weigh basis functions may performed by re-optimizing the DPD coefficients x according to:

${E_{V}(x)} = {{\frac{\rho}{n}{x}^{2}} + {\frac{1}{m}{\sum\limits_{i = 1}^{m}{{e_{i} - {\sum\limits_{k = 1}^{n}{x_{k}{B_{k}\left( {q_{i},p_{i}} \right)}}}}}^{2}}}}$ where m is the i-th column of the database matrix (i.e., the V matrix), and x_(k) (k=1, . . . , n) are the elements of a complex n-dimensionalized vector of the compensator (which, when adapted using the approaches described herein, acts as a universal compensator), with x∈

^(n).

To illustrate the effect of using a variety-based approach to select well distributed sample data points to perform DPD coefficient optimization processing, consider FIG. 6A showing graphs 610 and 612 depicting the simulation results for an optimization process performed without ensuring a good variety of samples. In this simulation, 100 samples of limited variety input/output data with argument values randomly spread over [0, 1] were generated. The distribution of the samples is shown in graph 610 of FIG. 6A. An optimization process to fit polynomials of degree 3 was then performed. The graph 612 shows the fitting error (mismatch) resulting from the use of the limited variety sample points. As expected, large fitting error, ranging over the [−1, 1] resulted from the fitting procedure applied to the limited variety samples.

In contrast, FIG. 6B shows simulation results when the optimization process used samples that were subject to a variety-based distribution approach that promoted sample diversity. In this simulation, new data was continually used to update a variety database matrix with 100 records (i.e., 100 columns). The first 1000 samples (with x and y arguments) that were fed and processed by the sample selection and replacement process (based on an approach similar to that described above) had a range of x∈[−1, 0], y ∈[−1, 0], then 1000 samples with x∈[−1, 0], y ∈[0, 1] were processed, and finally, 1000 samples with x∈[0, 1], y ∈[−1, 0] were processed. As shown in graph 620 of FIG. 6B, the coverage process (i.e., sample selection/replacement process) maintains a reasonably uniform coverage (and thus has high diversity/variety level) over the whole range [−1, 1]. Consequently, and as shown in graph 622, the fitting mismatch resulting from performing an optimization process (similar to that applied in relation to the simulation of FIG. 6A) is much improved over the resultant fitting mismatch of the simulation results illustrated in FIG. 6A.

With reference next to FIG. 7, a flowchart of an example procedure 700 for digital predistortion (particularly, to perform DPD adaptation using a variety-based approach) is shown. The procedure 700 includes collecting 710 data for a wireless device comprising a transmit chain (e.g., similar to the transmission section F 120 depicted in FIG. 1) that includes at least one power amplifier that produces output with non-linear distortions, with the data including data records that each includes input signal data (e.g., the signal u[t]) and corresponding output signal data (v[t] and/or y[t] signals shown in FIG. 1) for the transmit chain, and further being associated with a set of operating conditions values (shown as the vector p[t] in FIGS. 1 and 7). The each data record is associated with a 7variety measure representative of an extent of distinguishability of the respective data record from other collected data records. Examples of the operating conditions for the each data record may include one or more of, for example, channel information associated with a corresponding channel used by the wireless device to process the respective data record, a corresponding temperature, a corresponding average power, a corresponding power supply voltage (V_(DD)) associated with the at least one power amplifier of the transmit chain, a corresponding Voltage Standing Wave Ratio (VSWR) associated with the wireless device, and/or age of the each data record.

As noted, in some implementations, the predistorter may be realized in accordance with the first technique approach discussed above. In that approach, adaptation of the predistorter may be implemented based on samples that are collected in narrow subbands, rather than samples collected under wideband conditions (e.g., in which the bandwidth of the device's wireless transceiver obtaining the data is substantially equal to the spectral range for which the predistorter is to be implemented). Thus, in such embodiments, collecting the data for the wireless device may include collecting the data in one or more spectral channels, the one or more spectral channels located within a wide band spectral range, with each of the one or more spectral channels having a respective bandwidth smaller than a range bandwidth of the spectral range.

In some examples, each of the data records may include one or more samples obtained by the wireless device under a particular set of operating conditions. The each of the data records may be maintained in a separate column of a records matrix data structure.

With continued reference to FIG. 7, the procedure 700 further includes selecting 720 at least some of the data records according to a variety criterion computed based on the variety measure for the each data record. Selecting the at least some of the data records according to the variety criterion computed based on the variety measure for the each data record may include selecting data records associated with respective variety measures indicative of high distinguishability in order to increase diversity of the data records used to compute the digital predistortion parameters.

The procedure 700 additionally includes computing 730, based on the selected at least some of the data records, digital predistortion parameters used to control a digital predistorter of the wireless device, and applying 740 the digital predistorter, configured with the computed digital predistortion parameters, to subsequent input signals. In some embodiments, applying the digital predistorter, configured with the digital predistortion parameters, may include weighing basis functions of the digital predistorter, with the basis functions including (in embodiments such as those discussed in relation to the adaptive, expanded, basis functions that are formed from a base set of basis functions) complex terms applied to values of input samples provided in the data records, and with at least some of the basis functions further comprising real terms applied to values of sets of operating conditions associated with the data records.

In some implementations, the procedure 700 may also include computing, for the each data record, the variety measure based on distance metric representative of distances between data elements of a particular data record, to respective data elements for at least some of the collected data records. The computed variety measure for the each data record may be maintained in a distance vector data structure (such as the R vector discussed above) separate from a data structure to store the collected data records.

As discussed, under the variety approach, data samples used for the adaptation procedure are occasionally replaced with new sample data records that may be more relevant in that they promote diversity of the samples used in the adaptation process (in part because they are newer). Examples of a replacement procedure for data records includes occasionally collecting additional data records, storing a particular one of the collected additional data records in a records data structure in response to a determination that a variety measure computed for the particular one of the collected additional data records is more distinguishable, relative to previously stored data records available at the records data structure, than one of the previously stored data records, and re-computing, at regular or irregular intervals, the digital predistortion parameters based on the stored data samples available in the records data structure, including the stored particular one of the collected additional data records. In such embodiments, storing the particular one of the collected additional data records may include replacing in the records data structures the one of the previously stored data records with the particular one of the collected additional data records. The procedure may also include computing a corresponding variety measure for the particular one of the collected additional data records, including determining for each pair of records comprising the particular one of the collected additional data records and at least one of the previously stored data records a maximal absolute value of difference between the corresponding elements of the data records.

The approaches described above may be used in conjunction with the techniques described in PCT Application No. PCT/US2019/031714, filed on May 10, 2019, titled “Digital Compensation for a Non-Linear System,” which is incorporated herein by reference. For instance, the techniques described in that application may be used to implement the actuator (referred to as the pre-distorter in the incorporated application), and to adapt its parameters, and in particular to form the actuator to be responsive to an envelope signal or other signal related to power control of a power amplifier.

In some implementations, a computer accessible non-transitory storage medium includes a database (also referred to a “design structure” or “integrated circuit definition dataset”) representative of a system including, in part, some or all of the components of a digital predistortion, and/or DPD adaptation implementations and processes described herein. Generally speaking, a computer accessible storage medium may include any non-transitory storage media accessible by a computer during use to provide instructions and/or data to the computer. For example, a computer accessible storage medium may include storage media such as magnetic or optical disks and semiconductor memories. Generally, the database representative of the system may be a database or other data structure which can be read by a program and used, directly or indirectly, to fabricate the hardware comprising the system. For example, the database may be a behavioral-level description or register-transfer level (RTL) description of the hardware functionality in a high-level design language (HDL) such as Verilog or VHDL. The description may be read by a synthesis tool which may synthesize the description to produce a netlist comprising a list of gates from a synthesis library. The netlist comprises a set of gates which also represents the functionality of the hardware comprising the system. The netlist may then be placed and routed to produce a data set describing geometric shapes to be applied to masks. The masks may then be used in various semiconductor fabrication steps to produce a semiconductor circuit or circuits corresponding to the system. In other examples, the database may itself be the netlist (with or without the synthesis library) or the data set.

Unless defined otherwise, all technical and scientific terms used herein have the same meaning as commonly or conventionally understood. As used herein, the articles “a” and “an” refer to one or to more than one (i.e., to at least one) of the grammatical object of the article. By way of example, “an element” means one element or more than one element. “About” and/or “approximately” as used herein when referring to a measurable value such as an amount, a temporal duration, and the like, encompasses variations of ±20% or ±10%, ±5%, or +0.1% from the specified value, as such variations are appropriate in the context of the systems, devices, circuits, methods, and other implementations described herein. “Substantially” as used herein when referring to a measurable value such as an amount, a temporal duration, a physical attribute (such as frequency), and the like, also encompasses variations of ±20% or ±10%, ±5%, or +0.1% from the specified value, as such variations are appropriate in the context of the systems, devices, circuits, methods, and other implementations described herein.

As used herein, including in the claims, “or” as used in a list of items prefaced by “at least one of” or “one or more of” indicates a disjunctive list such that, for example, a list of “at least one of A, B, or C” means A or B or C or AB or AC or BC or ABC (i.e., A and B and C), or combinations with more than one feature (e.g., AA, AAB, ABBC, etc.). Also, as used herein, unless otherwise stated, a statement that a function or operation is “based on” an item or condition means that the function or operation is based on the stated item or condition and may be based on one or more items and/or conditions in addition to the stated item or condition.

Although particular embodiments have been disclosed herein in detail, this has been done by way of example for purposes of illustration only, and is not intended to be limit the scope of the invention, which is defined by the scope of the appended claims. Features of the disclosed embodiments can be combined, rearranged, etc., within the scope of the invention to produce more embodiments. Some other aspects, advantages, and modifications are considered to be within the scope of the claims provided below. The claims presented are representative of at least some of the embodiments and features disclosed herein. Other unclaimed embodiments and features are also contemplated. 

What is claimed is:
 1. A method for digital predistortion of signals provided to a radio frequency (RF) transmission path configured to transmit radio signals in a plurality of subbands within a spectral range, the method comprising: configuring a digital predistorter for predistorting signals comprising arbitrary spectral content within the spectral range, the configuring including: acquiring data samples representing operation of the RF transmission path to transmit radio signals in different subbands, each data sample including a digital input signal provided to the RF transmission path for transmission and a digital sensed signal representing a sensing of a radio transmission in response to the provided digital signal, the each digital input signal representing spectral content concentrated in a respective single corresponding subband within the spectral range, and updating parameters of the digital predistorter according to the acquired data samples to mitigate non-linear characteristics of the RF transmission path; receiving a further input signal for transmission by the RF transmission path, the further input signal representing spectral content in a particular subband within the spectral range; and using the configured predistorter to process the further input signal to yield a predistorted signal for providing to the RF transmission path.
 2. The method of claim 1, wherein the spectral range, W, comprises n RF bands of approximately equal-sized bandwidths.
 3. The method of claim 1, wherein acquiring the one or more data samples comprises: controllably adjusting a wireless device comprising the RF transmission path to operate in one or more of the different subbands within the spectral range; and at each of the one or more of the different subbands, measuring at least one set of samples.
 4. The method of claim 1, wherein acquiring the data samples comprises: measuring during runtime, subsequent to an initial computation of the parameters of the digital predistorter, additional sets of samples at additional subbands within the spectral range; and re-updating the parameters of the digital predistorter further based on at least some of the additional sets of samples.
 5. The method of claim 1, wherein the digital predistorter is implemented in a user equipment (UE), and wherein the plurality of subbands are successively assigned bands for transmission.
 6. The method of claim 1, wherein the updated parameters of the digital predistorter are configured to weigh non-linear functions used to predistort the further input signal, wherein the non-linear functions include complex terms that depend on a τ-dimensional complex vector q corresponding to complex input samples to the digital distorter, and real terms that depend on a d-dimensional real vector p representative of values of operating condition parameters associated with operation of the RF transmission path, wherein d is representative of a number of the operating condition parameters associated with the operation of the RF transmission path.
 7. The method of claim 6, wherein a total number of the non-linear functions used to digitally predistort input signals is proportional to the product of a number of terms, b, in a first set of non-linear functions, and d+1, such that the total number of non-linear functions is proportional to b·(d+1).
 8. The method of claim 1, wherein acquiring the data samples representing the operation of the RF transmission path in the different subbands further comprises: collecting for the one or more data samples corresponding operating condition parameters, wherein the corresponding operating condition parameters comprise one or more of: a corresponding temperature of the RF transmission path, a corresponding average power, a corresponding power supply voltage (V_(DD)) associated with at least one power amplifier of the RF transmission path, a corresponding Voltage Standing Wave Ratio (VSWR) associated with the RF transmission path, age of the collected each of the sets of samples, or corresponding channel information in which the data samples are collected.
 9. The method of claim 1, wherein using the configured predistorter to process the further input signal comprises: predistorting the further input signal according to: ${{v\lbrack t\rbrack} = {{u\lbrack t\rbrack} + {\sum\limits_{k = 1}^{n}{x_{k}{B_{k}\left( {{q_{u}\lbrack t\rbrack},{p\lbrack t\rbrack}} \right)}}}}},{{{where}\mspace{14mu}{q_{u}\lbrack t\rbrack}} = \begin{Bmatrix} {u\lbrack t + l - 1\rbrack} \\ \vdots \\ {u\left\lbrack {t + l - \tau} \right\rbrack} \end{Bmatrix}}$ where v[t] is an output of the digital predistorter at time t, u[t] is an input sample to the digital predistorter at the time t, B_(k) are the non-linear functions for k=1 to n, q_(u)[t] are stacks of input samples, p[t] are d-dimensional real vectors of operating condition parameters at the time t, and x_(k) are the computed digital predistortion parameters to weigh the non-linear functions.
 10. The method of claim 9, wherein updating the parameters of the digital predistorter comprises: deriving, according to an optimization process that is based, at least in part, on the acquired data samples measured in the different subbands, the digital predistortion parameters, x_(k), to weigh the non-linear functions B_(k).
 11. The method of claim 10, wherein deriving according to the optimization process the parameters, x_(k), of the digital predistorter comprises: adaptively adjusting the digital predistortion parameters, x_(k), through minimization of a regularized mean square error function E(x) according to: ${{E(x)} = \left. {{\frac{\rho}{n}{x}^{2}} + {\frac{1}{N}\sum\limits_{t \in T}}} \middle| {{v\lbrack t\rbrack} - {y\lbrack t\rbrack} - {\sum\limits_{k = 1}^{n}{x_{k}{B_{k}\left( {{q_{u}\lbrack t\rbrack},{p\lbrack t\rbrack}} \right)}}}} \right|^{2}},$ where ρ>0 is a regularization factor, v[t] is an output of the digital predistorter at the time t, and y[t] is the output of an observation receiver coupled to an output of the RF transmission path.
 12. The method of claim 11, wherein adaptively adjusting the digital predistortion parameters, x_(k), comprises: deriving correlation matrices G and L, when n<<100, with n being representative of the number of non-linear functions, according to: ${G = {\frac{1}{N}{\sum\limits_{t \in T}{{B\lbrack t\rbrack}{B\lbrack t\rbrack}^{\prime}}}}},{L = {\frac{1}{N}{\sum\limits_{t \in T}{\left( {{v\lbrack t\rbrack} - {y\lbrack t\rbrack}^{\prime}} \right){B\lbrack t\rbrack}}}}},{{{{where}\mspace{14mu}{B\lbrack t\rbrack}} = \begin{bmatrix} {B_{1}\left( {{q_{y}\lbrack t\rbrack},{p\lbrack t\rbrack}} \right)} \\ \vdots \\ {B_{n}\left( {{q_{y}\lbrack t\rbrack},{p\lbrack t\rbrack}} \right)} \end{bmatrix}};}$ and deriving the digital predistortion parameters, x_(k), according to: $\overset{\_}{x} = {\left( {{\left( \frac{p}{n} \right)I} + G} \right)^{- 1}{L.}}$
 13. The method of claim 11, wherein adaptively adjusting the digital predistortion parameters, x_(k), comprises: approximating the digital predistortion parameters, x_(k), according to a stochastic gradient process.
 14. The method of any of claim 1, wherein acquiring the data samples representing operation of the RF transmission path to transmit radio signals in the different subbands comprises: selecting the data samples from a plurality of data samples according to a variety criterion computed based on the variety measure for the each data record, wherein the plurality of records are arranged in a variety database structure in which each record corresponds to one of the data samples.
 15. A linearization system to perform digital predistortion of signals provided to a radio frequency (RF) transmission path configured to transmit radio signals in a plurality of subbands within a spectral range, the system comprising: the RF transmission path comprising at least one power amplifier that produces output with non-linear distortions; and a digital predistorter comprising an adaptation section and a compensator, the digital predistorter configured to: configure the compensator to predistort signals comprising arbitrary spectral content within the spectral range, including: acquire data samples representing operation of the RF transmission path to transmit radio signals in different subbands, each data sample including a digital input signal provided to the RF transmission path for transmission and a digital sensed signal representing a sensing of a radio transmission in response to the provided digital signal, the each digital input signal representing spectral content concentrated in a respective single corresponding subband within the spectral range, and update parameters of the digital predistorter according to the acquired data samples to mitigate non-linear characteristics of the RF transmission path; receive a further input signal for transmission by the RF transmission path, the further input signal representing spectral content in a particular subband within the spectral range; and use the configured compensator to process the further input signal to yield a predistorted signal for providing to the RF transmission path.
 16. The system of claim 15, wherein the digital predistorter configured to acquire the data samples is configured to: controllably adjust a wireless device comprising the RF transmission path to operate in one or more of the different subbands within the spectral range; and at each of the one or more of the different subbands, measure at least one set of samples.
 17. The system of claim 15, wherein the digital predistorter configured to acquire the data samples is configured to: measure during runtime, subsequent to an initial computation of the parameters of the digital predistorter, additional sets of samples at additional subbands within the spectral range; and re-update the parameters of the digital predistorter further based on at least some of the additional sets of samples.
 18. The system of claim 15, wherein the digital predistorter configured to use the compensator is configured to: derive from a first set of non-linear functions an expanded set of non-linear functions of the digital predistorter, the expanded set of non-linear functions including (d+1) groups of non-linear functions, with d being representative of a number of operating condition parameters, wherein a first group of the expanded set of non-linear functions includes the first set of non-linear functions without any modification, and wherein each of the remaining d groups of the expanded set of non-linear functions is determined as a transformation of the first set of non-linear functions according to a function of at least one of the operating condition parameters; and weigh the expanded set of non-linear functions with the updated digital predistortion parameters, to process the further input signal.
 19. The system of claim 15, wherein the digital predistorter configured to acquire the data samples representing operation of the RF transmission path to transmit radio signals in the different subbands is configured to: select the data samples from a plurality of data samples according to a variety criterion computed based on the variety measure for the each data record, wherein the plurality of records are arranged in a variety database structure in which each record corresponds to one of the one or more data samples.
 20. A non-transitory machine-readable medium storing a design structure comprising elements that, when processed in a computer-aided design system, generate a machine-executable representation of a linearization system that is used to fabricate hardware comprising: a radio frequency (RF) transmission path circuit comprising at least one power amplifier that produces output with non-linear distortions; and a digital predistorter circuit comprising an adaptation section and a compensator, the digital predistorter configured to: configure the compensator to predistort signals comprising arbitrary spectral content within the spectral range, including: acquire data samples representing operation of the RF transmission path circuit to transmit radio signals in different subbands, each data sample including a digital input signal provided to the RF transmission path circuit for transmission and a digital sensed signal representing a sensing of a radio transmission in response to the provided digital signal, the each digital input signal representing spectral content concentrated in a respective single corresponding subband within the spectral range, and update parameters of the digital predistorter according to the acquired data samples to mitigate non-linear characteristics of the RF transmission path circuit; receive a further input signal for transmission by the RF transmission path circuit, the further input signal representing spectral content in a particular subband within the spectral range; and use the configured compensator to process the further input signal to yield a predistorted signal for providing to the RF transmission path circuit. 